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by "plane waves" in terms of the corresponding radial coordinates. We also present an alternative 
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relations generalize those derived earlier for the multifractal spectrum of the leading operators. 
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I. INTRODUCTION 

The phenomenon of Anderson localization of a quan- 
tum particle or a classical wave in a random environ- 
ment is one of the central discoveries made by condensed 
matter physics in the second half of the last century.EJ 
Although more than fifty years have passed since An- 
derson's original papcs, Anderson localization remains a 
vibrant research field.BI One of its central research direc- 
tions is the physics of Anderson transitions,™ including 
metal-insulator transitions and transitions of quantum- 
Hall type (i.e. between different phases of topological in- 
sulators). While such transitions are conventionally ob- 
served in electronic conductor and semiconductor struc- 
tures, there is also a considerable number of other experi- 
mental realizations actively studied in recnat and currertL 
works. Those include localization of light,Q cold atoms,El 
ultrasoundjj and optically driven atomic systems.Q On 
the theory side, the field received a strong boost through 
the discovery of unconventional symmetry classes and the 
development of ar-fjomplete symmetry classification of dis- 
ordered systems. clErfcj The unconventional classes emerge 
due to additional particle- hole and/or chiral symmetries 
that are, in particular, characteristic for models of dis- 
ordered superconductors and disordered Dirac fermions 
(e.g. in graphene). In total one has 10 symmetry classes, 
including three standard (Wigner-Dyson) classes, three 
chiral, and four Bogoliubov-de Gennes ("superconduct- 
ing") classes. This multitude is further supplemented by 
the possibility for the underlying field theories to have 
a non-trivial topology (9 and Wess-Zumino terms), lead- 



ing to a rich "zoo" of Andersoi|t4;ransition critical points. 
The recent advent of graphenelHl and of topological insu- 
lators and superconductorala reinforced the experimental 
relevance of these theoretical concepts. 

In analogy with more conventional second-order phase 
transitions, Anderson transitions fall into different uni- 
versality classes according to the spatial dimension, sym- 
metry, and topology. In each of the universality classes, 
the behavior of physical observables near the transition 
is characterized by critical exponents determined by the 
scaling dimensions of the corresponding operators. 

A remarkable property of Anderson transitions is that 
the critical wave functions are multifractal due to their 
strong fluctuations. Specifically, the wave-function mo- 
ments show anomalous multifractal scaling with respect 
to the system size L, 

i'^(|V(r)P«)(xL-^^ T, ^ d{q - 1) + A,, (1.1) 

where d is the spatial dimension, (. . .) denotes the oper- 
ation of disorder averaging and are anomalous multi- 
fractal exponents that distinguish the critical point from 
a simple metallic phase, where = 0. Closely related 
is the scaling of moments of the local density of states 
(LDOS) i^(r), 

(l^-?) CX Xg=Ag+qX^, (1.2) 

where Xi, = x\ controls the scaling of the average LDOS, 
{v) oc L~^" . Multifractality implies the presence of 
infinitely many relevant (in the renormalization-group 
(RG) sense) operators at the Anderson-transition critical 
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point. First steps towards the experimental deteriniiia^. 
tion of multifractal spectra have been made recentlyJ3l3'E£l 

Let us emphasize that when we speak about a q-th 
moment, we neither require that q is an integer nor that 
it is positive. Throughout the paper, the term "moment" 
is understood in this broad sense. 

In Refs. a symmetry relation for the LDOS dis- 

tribution function (and thus, for the LDOS moments) in 
the Wigner-Dyson symmetry classes was derived: 



(1.3) 



with = 1. Equation (L^) is obtained in the framework 
of the non-linear cr-model and is fully general otherwise, 
i.e., it is equally applicable to metallic, localized, and 



critical systems. An important consequence of Eq. (1.3) 
is an exact symmetry|-^elation for Anderson-transition 
multifractal exponent; 



-q- 



(1.4) 



While (T-models in general are app roximations to partic- 
ular microscopic systems, Eq. (1.4) is exact in view of the 
universality of critiftal behavior. 

In a recent paper ,113 the three of us and A. W. W. Lud- 
wig uncovered t he g ro up- theoretical origin of the sym- 
metry relations (|l.3|), (|l.4D. Specifically, we showed that 
these relations are manifestations of a Weyl symmetry 
group acting on the cr-modcl manifold. This approach 
was further used to generalize these relations to the un- 
conventional (Bogoliubov-de Gennes) classes CI and C, 
with = 2 and = 3, respectively. 

The ope rators representing the averaged LDOS mo- 
ments (1.2) by no means exhaust the composite operators 
characterizing LDOS (or wave-function) correlations in a 
disordered system. They are distinguished in that they 
arc the dominant (or most relevant) operators for each 
but they only represent "the tip of the iceberg" of a much 
larger family of gradientless composite operators. Often, 
the subleading operators are also very important phys- 
ically. An obvious example is the two-point correlation 
function 

AXri,r2) = |^^(ri)V'|(r2)| 

- Vo(ri)#(r2)C(r2)V';(ri), (1.5) 

which enters in the Hartree-Fock matrix element of a 
two-body interaction, 

M„0 = j dridr2Ka,p{Ti,Y2)U{Ti-V2). (1.6) 

Questions about the scaling of the disorder-averaged 
function i^Q,^ (ri, r2), its moments, and the correla- 
tions of such objects, arise naturally when one studies, 
e.g., the interaction-inducjcui.dephasing at the Anderson- 
transition critical point .Eto 

The goals and the results of this paper are threefold: 

1. We develop a systematic and complete classifica- 



tion of gradientless composite operators in the su- 
persymmetric non-linear cr-models of Anderson lo- 
calization. Our [approach here-difiers from that of 
Hof and WegnerEJ and Wegneifi3c3 in two respects. 
Firstly, we work directly with the supersymmetric 
(SUSY) theories rather than with their compact 
replica versions as in Refs. plj-p^. Secondly, we 
employ (a superization of) the Iwasawa decompo- 
sition and the Harish- Chandra isomorphism, which 
allow us to explicitly construct "radial plane waves" 
that are eigenfunctions of the Laplace-Casimir op- 
erators of the CT-model symmetry group, for arbi- 
trary (also non-integer, negative, and even com- 
plex) values of a set of parameters [generalizing 
the order q of the moment in Eq. (|l . l|) , (|l.2| )]. We 
also develop a more basic construction of scaling 
operators as highest-weight vectors (and explain 
the link with the Iwasawa-decomposition formal- 
ism). 

2. We establish a connection between these compos- 
ite operators and the physical observables of LDOS 
and wave-function correlators, as well as with some 
transport observables. 

3. Furthermore, the Iwasawa-decomposition formal- 
ism allows us to exploit a certain Weyl-group in- 
variance and deduce a large number of relations be- 
tween the scaling dimensions of various composite 
operators at criticality. These symmetry relations 
generalize Eq. (1.4) obtained earlier for the most 
relevant operators (LDOS moments). 

It should be e mphasized that we do not attempt to gen- 
eralize Eq. (1.3), which is also valid away from criticality, 
but rather focus on Anderson-transition critical points. 
The reason is as follows. The derivation of Eq. ( |l.3| ) 
in Ref. jl7| was based on a (super-)generalization of a 
theorem due to Harish-Chandra. We are not able to fur- 
ther generalize this theorem to the non -min imal a-models 
needed for the generalization of Eq. (1.3) to subleading 
operators. For this reason, we use a weaker version of 
the Weyl-invariance argument which is applicable only 
at criticality. This argument is sufficient to get exact 
relations between the critical exponents. 

In the main part of the paper we focus on the unitary 
Wigner-Dyson class A. Generalizations to other symme- 
try classes, as well as some of their peculiarities, are dis- 
cussed at the end of the paper. 

The structure of the paper is as follows. In Sec. || we 
briefly review Wegner's classification of composite oper- 
ators in replica cr-models by Young diagrams. In Sec. |l| 
we introduce the Iwasawa decomposition for supersym- 
metric cr-models of Anderson localization and, on its ba- 
sis, develop a classification of the composite operators. 
The correspondence between the replica and SUSY for- 
mulations is established in Secjvl for the case of the 
minimal SUSY model. Section 01s devoted to the con- 
nection between the physical observables (wave-function 
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correlation functions) and the cr-model composite opcr 
ators. This subject is further developed in Sees. VIA 



VI D| , where we identify observables that correspond to 



exact scaling operators and thus exhibit pure power scal- 
ing (without any ad mixtur e of subleading power-law con- 
tributions). In Sec. VIE we formulate a complete ver- 
sion (going beyond the minimal-SUSY model considered 
in Sec. IV) of the correspondence between the full set 
of operators of our SUSY classification and the physical 
observables (wave-function and LDOS correlation func- 
tions). An alternative and more basic approach to scaling 
operators via the notion of highest-weight vector is ex- 
plained in Sec. VI F, We also indicate how this approach 



is relat ed to the one based on the Iwasawa decomposition. 
In Sec. VII we employ the Weyl-group invariance and de- 
duce symmetry relations among the anomalous dimen- 
sions of various composite operators at criticality. The 
generalization of t hese r esults to other symmetry classes 
is discussed in Sec. VIII . In Sec. [X we analyze the impli- 
cations of our findings for transport observables defined 
within the Dorokhov-Mello-Pereyra-Kumar (DMPK) for- 
malism. In Sec. ^ we discuss peculiarities of symmetry 
classes whose cr- models possess additional 0(1) [classes 
D and Dili] or U(l) [classes BDI, CII, Dili] degrees of 
freedom. Section XI contains a summary of our results, 



as well as a discussion of open questions and directions 
for further research. 



II. REPLICA cr-MODELS AND 
WEGNER'S RESULTS 

The replica method leads to the reformulation of the 
localization problem as a theory of fields taking ¥al4es 
in a symmetric space G/K - a non-linear a-model.ac3 If 
one uses fermionic replicas, the resulting a-model target 
spaces are compact, and for the Wigner-Dyson unitary 
class (a.k.a. class A) they are of the type G/K with G = 
U(toi -I- TO2) and K = U(mi) x U(to2). Bosonic replicas 
lead to the non-compact counterpart G' / K where G' = 
U(mi,m2). The total number of replicas m = mi + m2 
is taken to zero at the end of any calculation, but at 
intermediate stages it has to be sufficiently large in order 
for the (T-model to describe high enough moments of the 
observables of interest. The cr-model field Q is a matrix, 
Q = gAg~^, where A = diag(lmi , — Ima) g £ G. 
Since Q does not change when g is multiplied on the 
right {g — >■ gk) by any element k G K, the set of matrices 
Q realizes the symmetric space G/K. Clearly, Q satisfies 
the constraint ^ 1. Roughly speaking, one may think 
of the symmetric space G/K as a "generalized sphere". 

The action functional of the cr-model has the following 
structure: 

S[Q] = ^ j d'^r Tr(VQ)2 + h J d'^r TvQA. (2.1) 

Here Q(r) is the Q- matrix field depending on the spatial 
coordinates r. The parameter l/lGnt in front of the first 



term is 1/8 times the conductivity (in natural units). In 
the renormalization-group (RG) framework, t serves as a 
running coupling constant of the theory. While the first 
term is invariant under conjugation Q{r) .9oQ(r).9(7^ 
of the Q-matrix field by any (spatially uniform) element 
go G G, the second term causes a reduction of the symme- 
try from G to K, i.e. only conjugation Q(r) kQQ{r)kQ^ 
hy ko £ K leaves the full action invariant. The second 
term provides an infrared regularization of the theory 
in infinite volume; physically, h is proportional to the 
frequency. When studying scaling properties, it is usu- 
ally convenient to work at an imaginary frequency, which 
gives a non-zero width to the energy levels. If the physical 
system has a spatial boundary with coupling to metallic 
leads, then boundary terms arise which are X-invariant 
like the second term in (2.1). 



Quite generally, physical observables are represented 
by composite operators of the corresponding field the- 
ory. For the case of the compact cr-model resulting from 
fermionic replicas, a classification of composite operators 
withoutj— spatial derivatiMcs, was developed by Hof and 
Wegnercll and WegncrBSEj It goes roughly as follows. 
The composite operators were constructed as polynomi- 
als in the matrix elements of Q, 



P = 



Qh ,12 ' ' ' Qi 



(2.2) 



Such polynomials transform as tensors under the action 
{Q ~^ 9Q9~^) of the group G — U(mi + m^)- They 
decompose into polynomials that transform irreducibly 
under G, and composite operators (or polynomials in 
Q) belonging to different irreducible representations of 
the symmetry group G do not mix under the RG flow. 
The renormalization within each irreducible representa- 
tion is characterized by a single renormalization con- 
stant. Therefore, fixing any irreducible representation 
it is sufficient to focus on operators in a suitable one- 
dimensional subspace. In view of the iC-symmetry of 
the action, a natural choice of subspace is given by K- 
invariant operators, i.e. those polynomials that satisBz. 
P[Q) = P{koQk^^) for aU ko € K. It can be shownEl 
that each irreducible representation occurring in ( |2.2| ) 
contains exactly one such operator. We may therefore 
restrict our attention to A'-invariant operators. By their 
A'-invariancc, such operators can be represented as linear 
combinations of operators of the form 



Px = Tr(AQ)''^i • • • Tr(AQ)^ 



(2.3) 



where £ ~ minjmi, TO2}, and A = {fci, . . . , ki} is a parti- 
tion k = ki + . . . + ki such that ki > . . . > ki > 0. In 
particular, for fc = 1 we have one such operator, {1}, for 
k ^ 2 < £ two operators, {2} and {1, 1}, for fc = 3 < ^ 
three operators {3}, {2, 1}, and {1, 1, 1}, for fc = 4 < £ 
five operators {4}, {3, 1}, {2, 2}, {2, 1, 1}, and {1, 1, 1, 1}, 
and so on.l2J As described below, operators of order fc cor- 
respond to observables of order fc in the LDOS or. equiv- 
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alently, of order 2k in the wave-function amplitudes. 

It turns out that the counting of partitions yields the 
number of different irreducible representations that oc- 
cur for each order k of the operator. More precisely, 
there is a one-to-one correspondence between the irre- 
ducible representations of G ~ U(mi -I- 7712) which occur 
in ( p..2\ ) and the set of irreducible representations of U(£), 
£ = min(?Tii, 7712), as given by partitions A = (fci, . . . , ki). 
We may also think of the partition A Young dia- 
gram for U(£). (Please note that Young diagrams and 
the corresponding partitions are commonly denoted by 
using parentheses as opposed to the curly braces of the 
above discussion. An introduction to Young diagrams 
and their use in our context is given in Appendix |AJ.) 

The claimed relation with the representation theory of 
U(^) becomes plausible if one uses the Cartan decompo- 
sition G ~ KAK, by which each element of G is repre- 
sented as g = kak', where k,k' G K, a G A, and A ~ 
U(l)^ is a maximal abelian subgroup of G with Lie alge- 
bra contained in the tangent space of G/K at the origin. 
In this decomposition one has Q = kaAa~^k~^ . A K- 
invariant operator P satisfies P{Q) = P{kaAa^^k^^) = 
P{aAa~^). In other words, P depends only on a set of 
£ "ii'-radial" coordinates for a G A ~ U(l)^ - this is 
ultimately responsible for the one-to-one correspondence 
with the irreducible representations of U(£). 

The X-invariant operators associated with irreducible 
representations are known as zonal spherical functions. 
For the case of G/K = U(2)/U(l) x U(l) = 5*2, which is 
the usual two-sphere, they are just the Legendre polyno- 
mials, i.e. the usual spherical harmonics with magnetic 
quantum number zero; see also Appendix Please note 
that here and throughout the paper we use the conven- 
tion that the symbol for the direct product takes prece- 
dence over the symbol for the quotient operation. Thus 

G/Ki X A'2 = G/{Ki XK2). 

From the work of Harish-Chandrall it is known that the 
zonal spherical functions have a very simple form when 
expressed by A'^-radial coordinates that originate from the 
Iwasawa decomposition G = NAK; see Sec. [II below. 



This will make it possible to connect Wegner's classifica- 
tion of composite operators with our SUSY classification, 
where we use the Iw^asawa decomposition. 

Hof and Wegncrnil calculated the anomalous dimen- 
sions of the polynomial composite operators ( |2.2| ) for a- 
models on the target spaces G{mi + 1112) / G{mi) x G{m2) 
for G = O, U, and Sp (whose replica limits correspond to 
the Anderson localization problem in the Wigner-Dyson 
classes A, AI, and All, respectivdjdin 2 + e dimensions 
up to three-loop order. WegnerEScJ extended this cal- 
culation up to four-loop order. The results of Wegner 
for the anomalous dimensions are summarized in the C- 
f unction for each composite operator: 

Ut) = a2{X)p{t) + C(3)c3(A)i4 + 0{t'), (2.4) 

where t, serving as a small parameter of the expansion, is 



|A| 


A 


0-2 


C3 


1 


(1) 








2 


(2) 


2 


6 




(1,1) 


-2 


6 


3 


(3) 


6 


54 




(2,1) 










(1,1,1) 


-6 


54 


4 


(4) 


12 


216 




(3,1) 


4 


24 




(2,2) 










(2,1,1) 


-4 


24 




(1,1,1,1) 


-12 


216 


5 


(5) 


20 


600 




(4,1) 


10 


150 




(3,2) 


4 


24 




(3,1,1) 










(2,2,1) 


-4 


24 




(2,1,1,1) 


-10 


150 




(1,1,1,1,1) 


-20 


600 



TABLE I. CoefScients a2 and C3 of the (^-function for class A 
in the replica limit. Results for composite operators charac- 
terized by Young diagrams up to size |A| = fc = 5 are shown. 



the renormalized coupling constant of the cr-model. The 
coefficients 02 and C3 depend on the operator Pa (defined 
by the Young diagram A) and on the type of model (O, U, 
or Sp, as well as r?7,i and 7712). The function p(t) depends 
on the model only and not on A. The coefficient 02 hap- 
pens to be the quadratic Casimir eigenvalue ass ociated to 
the representation with Young diagram A (A6). For the 
case of unitary symmetry (class A), on which we focus, 
the coefficients satisfy the following symmetry relations: 



a2(A, m) -a2(A, -777) 
C3(A, 777) C3(A,-r77), 



(2.5) 
(2.6) 



where A is the Young diagram conjugate to A, i.e., A 
is obtained by reflection of A with respect to the main 
diagonal. By using the results of Table 2 from ReL , 
complementing them with these symmetry relations, and 
taking the replica limit 7771 = 7772 = 0, we can obtain the 
values of the coefficients 02 and C3 for all polynomial 
composite operators up to order fc = 5. These values are 
presented in Table ^. The function p{t) is given for this 
model (unitary case, replica limit) by 



3 .3 



p{t) ^t+^t 



(2.7) 



A note on conventions and nomenclature is in order 
here. The way we draw Young diagrams (see Appendix 
^) is the standard way. Thus the horizontal direction 
corresponds to symmetrization and the vertical one to an- 
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tisymmetrization. In Wegner's approach fcrmionic repli- 
cas are used, hence his natural observables are antisym- 
metrized products of wave functions, whereas the descrip- 
tion of symmetrized products (like LDOS moments) re- 
quires the symmetry group to be enlarged. Wegner uses 
a different convention for labeling the invariant scaling 
operators, employing the Young diagrams conjugate to 
the usual ones used here. Thus, for example, the LDOS 
moment {v'^) corresponds in our convention to the Young 
diagram (g), while it is labeled by (f^) in Wegner's works. 
This has to be kept in mind when comparing our Table 
I with Table 2 of Rcf . . Of course, if one uses bosonic 
replicas, the situation is reversed: the natural objects 
then are symmetrized products and the roles of the hor- 
izontal and vertical directions get switched. 

While the works pl|-|2^ signified a very important 
advance in the theory of critical phenomena described 
by non-linear cr-models, the classification of gradicntlcss 
composite operators developed there is complete only for 
compact models. This can be understood already by in- 
specting the simple example of U(2)/U(l) x U(l) = S'^ 
(two-sphere), which is the target space of the conven- 
tional 0(3) non-linear cr-model. As was mentioned above, 
the corresponding isT-invariant composite operators are 
the usual spherical harmonics Yjo with / = 0,1,..., 
which arc Legcndre polynomials in cos 9. (The polar 
angle 9 parametrizes the abelian group A, which is one- 
dimensional in this case.) It is well known that the spher- 
ical harmonics indeed form a complete system on the 
sphere. The angular momentum I plays the role of the 
size k = |A| of the Young diagram. The situation changes, 
however, when we pass to the non-compact counterpart, 
U(l, 1)/U(1) X U(l), which is a hyperboloid . The dif- 
ference is that now the polar direction (parametrized by 
the coordinate 9) becomes non-compact. For this reason, 
nothing forces the angular momentum I to be quantized. 
Indeed, the spherical functions on a hyperboloid arc 
characterized by a continuous parameter (determining 
the order of an associated Legendrc function) which takes 
the role of the discrete angular momentum on the sphere 

. See Appendix ^for more details. 

The above simple example reflects the general situ- 
ation: for theories defined on non-compact symmetric 
spaces the polynomial composite operators by no means 
exhaust the set of all composite operators. In the field 
theory of Anderson localization, we are thus facing the 
following conundrum: two theories, a compact and a 
non-compact one [U(mi -I- m2)/U(mi) x U(m2) resp. 
U(toi, TO2)/U(toi) X U(m2) for class A], which should 
describe in the replica limit mi = m2 = the same 
Anderson localization problem, have essentially different 
operator content. This is a manifestation of the fact that 
the replica trick has a very tricky character indeed. In 
this paper we resolve this ambiguity by using an alterna- 
tive, well-defined approach based on supersymmetry. 



III. SUSY a-MODELS: IWASAWA 
DECOMPOSITION AND CLASSIFICATION OF 
COMPOSITE OPERATORS 

In the SUSY formalism the cr-model target space is the 
coset space 

G/K = U(n,n|2n)/U(n|n) x \]{n\n). (3.1) 

This manifold combines compact and non-compact fea- 
tures "dressed" by anticommuting (Grassmann) vari- 
ables. Its base manifold Mq x Mi is a product of 
non-compact and compact symmetric spaces: Mq = 
U(ri,n)/U(n) x U(n) and A/j U(2n)/U(n) x U(n). 

The acti on fu nctional of the SUSY theory still has the 
same form (2.1), except that the trace Tr is now replaced 
by the supertrace STr. It is often useful to consider a 
lattice version of the model (i.e. with discrete rather than 
continuous spatial coordinates); our analysis based solely 
on symmetry considerations remains valid in this case 
as well. Furthermore, it also applies to models with a 
topological term (e.g., for quantum Hall systems in 2D). 

The size parameter n of the supergroups involved needs 
to be sufficiently large in order for the model to contain 
the observables of interest; this will be discussed in detail 
in Sec. The minimal variant of the model with n ~ 1 
can accommodate arbitraiai, moments {v"^) of the local 
density of states (LDOS) 1^,1111 but is in general insufficient 
to give more complex observables, e.g. moments of the 
Hartree-Fock matrix element (1^). We will first dospribe 
the construction of operators for the n = 1 modelllll and 
then the generalization for arbitrary n. 

Our approach is based oru^he Iwasawa decomposition 
for symmetric superspaces,E3'E3 generalizing the corre- 
sponding construction for non-compact classical symmet- 
ric spacesBa The Iwasawa decomposition factorizes G as 
G = NAK, where A is (as above) a maximal abelian 
subgroup for G/K, and A'^ is a nilpotent group defined 
as follows. One considers the adjoint action (i.e. the ac- 
tion by the commutator) of elements of the Lie algebra a 
of A on the Lie algebra g of G. Since a is abelian, all its 
elements can be diagonalized simultaneously. The cor- 
responding eigenvectors in the adjoint representation are 
called root vectors, and the eigenvalues are called roots. 
Viewed as linear functions on a, roots lie in the space 
a* dual to o. A system of positive roots is defined by 
choosing some hyperplane through the origin of a* which 
divides o* in two halves, and then defining one of these 
halves as positive. All roots that lie on the positive side 
of the hyperplane are considered as positive. The nilpo- 
tent Lie algebra n is generated by the set of root vectors 
associated with positive roots; its exponentiation yields 
the group N. The Iwasawa decomposition G — NAK 
represents any clement g e G in the form g = nak, with 
n & N , a & A, and k € K. This factorization is unique 
once the system of positive roots is fixed. 

An explanation is in order here. The Iwasawa decom- 
position G — NAK is defined as such only for the case 
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of a non-compact group G with maximal compact sub- 
group K. Now the latter condition appears to exclude the 
symmetric spaces G/K that arise in the SUSY context, as 
their subgroups K fail to be maximal compact in general. 
This apparent difficulty, however, can be circumvented by 
a process of analytic continuation^] Indeed, the classical 
Iwasawa decomposition G = NAK determines a triple 
of functions n : G ^ N, a : G ^ A, k : G ^ K hy 
the uniqueness of the factorization g = n(g)a{g)k{g) . In 
our SUSY context, where K is not maximal compact and 
the Iwasawa decomposition does not exist, the functions 
71(17), 0(5), and k[g) still exist, but they do as functions 
on G with values in the complexified groups Nc^ Ac, and 
Kic, respectively. In particular, the Iwasawa decomposi- 
tion gives us a multi- valued function In a which assigns 
to every group clement g G G an clement Ina(g) of the 
(complexification of the) abclian Lie algebra a. 

Note that for uq E N, ko G K one has a{nQgko) = 
a{g) by construction. Thus one gets a function a(Q) on 
G/K by defining d{gAg~^) = a{g). This function is A^- 
radial, i.e., it depends only on the "radial" factor A in the 
paramctrization G/K ~ NA and is constant along the 
nilpotent group N: alnoQriQ^) = a{Q). Its multi-valued 
logarithm In a{Q) will play some role in what follows. 

In the case n = 1, which was considered in Rcf. [ pT[ , 
the space a* is two-dimensional, and we denote its ba- 
sis of linear coordinate functions by x and y, with x 
corresponding to the boson-boson and y to the fermion- 
fermion sector of the theory. In terms of this basis we 
can choose the positive roots to be 

2x(l), 2iy(l), x + iy{^2), x-iy{^2), (3.2) 

where the multiplicities of the roots are shown in paren- 
theses; note that odd roots are counted with negative 
multiplicity. (A root is called even or odd depending 
on whether the corresponding eigcnspacc is in the even 
or odd part of the Lie superalgebra. Even root vectors 
are located within the boson-boson and fermion-fermion 
supermatrix blocks, whereas odd root vectors belong to 
the boson-fermion and fermion-boson blocks.) For this 
choice of positive root system, the Weyl co-vector p (or 
half the sum of the positive roots with multiplicities) is 



is constrained by 



-X + iy. 



(3.3) 



The crucial advantage of using the symmetric-space 
paramctrization generated by the Iwasawa decomposi- 
tion is that the iV-radial spherical functions (/j^ have the 
very simple form of exponentials (or "plane waves" ) , 

^ g(-l+M)x(lna(Q)) + (l+Mi)iy(lna(Q))^ ^g^^-j 

labeled by a weight vector fi — (Iqx + fJ-iiy in a*. The 
boson-boson component /io of the weight /i can be any 
complex number, while the fermion-fermion component 



yUi e {-1, -3, -5, . . .} 



(3.5) 



to ensure that g*(i+/^i)(ina(Q)) single-valued in spite of 
the presence of the logarithm. 

From here on we adopt a simplified notation where 
we use the same symbol x also for the composition of 
X with In a (and similar for y). Thus x may now have 
two different meanings: either its old meaning as a linear 
function on a, or the new one as an TV-radial function 
X o Ina on G/K. It should always be clear from the 
context which of the two functions x we mean. 

With this convention, Eq. ( |3.4D reads (^^ = e''+^ — 
g(-i+/.o):i+(i+Mi)*y. We will also use the notation 



1 



Mo 



1 



P = 



Ml 



e Z4 



(3.6) 



where Z+ means the set of non-negative integers. In this 



notation the exponential functions (3.4) take the form 



'Pp. 



2ipy 



(3.7) 



We mention in passing that the quantization of p is 
nothing but the familiar quantization of the angular 
momentum I for the well-known spherical functions on 

. Indeed, the "momentum" p is conjugate to the "ra- 
dial variable" y corresponding to the compact (fermion- 
fermion) sector. The absence of any quantization for q 
should also be clear from the discussion at the end of 
Sec. ||. In fact, q is conjugate to the radial variable x of 
the non-compact (boson-boson) sector, which is a hypcr- 
boloid H'^. 

By simple reasoning based on the observation that A 
normalizes N (i.e., for any oq G A and hq Cz N one has 
aQ^nottQ € N), each plane wave (/9^ is an eigenfunction 
of the Laplace-Beltrami operator ^and all other invari- 
ant differential operators on G/K!EB [The same conclu- 
sion follows from more general considerations based on 
highest-weight vectors (see Sec. VIF).] The eigenvalue 
of the Laplace-Beltrami operator is 



/io-M?=4g(g-l)-4p(p+l), 



(3.8) 



up to a constant factor. 

It should be stressed that the A''-radial spherical func- 
tions, which depend only on a in the Iwasawa decompo- 
sition g = nak, differ from the A'-radial spherical func- 
tions (depending only on a' in the Cartan decomposi- 
tion g — k'a'k" , see Sec. H), since for a given element 
Q = gA.g^^ or gK of G/K the radial elements a and a' of 
these two factorizations are different. However, a link be- 
tween the two types of radial spherical function can easily 
be established. Indeed, if ip{Q) is a spherical function, 
then for any element k £ K the transformed function 
ip{kgk~^) is still a spherical function from the same rep- 
resentation. Therefore, we can construct a A'-invariant 
spherical function by simply averaging Lp^{k~^Qk) 
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over K , 



up to a constant factor. 



JK 



(3.9) 



provided, of course, that the integral does not vanish. 

For n > 1 the space o* has dimension 2n. We label the 
linear coordinates with i ~ 1, . . . , n; following 

the notation above, the Xj and yj correspond to the non- 
compact and compact sectors, respectively. The positive 
root system can be chosen as follows: 



Xk (2), 



x-j + Xk (2), 



2x, (1), 



i(2/z - ym) (2), iiyi + ym) {2), 2iyi {l), 
Xj+iyi{~2), xj-iyi{^2), (3.10) 

where \ < j < k < n and 1 < m < I < n. As before, the 
multiplicities of the roots arc given in parentheses, and a 
negative multiplicity means that the corresponding root 
is odd, or fermionic. The half-sum of these roots (still 
weighted by multiplicities) now is 



P = ^CjXj + i'^biyi 



with 



= 1 - 2j, 



b, = 21- 1. 



(3.11) 



(3.12) 



The A^-radial spherical functions are constructed just 
like for n = 1. They are still "plane waves" (p^ = e^^^ 
but now the weight vector ^ has 2n components /i^ and 
n] , the latter of which take values 



III e {-fe(,-6,-2,-fe,-4,...} 
We will also write 

^l\ + h 



(3.13) 



Pi 



2 ' " 2 
In this notation our iV-radial spherical functions are 



(3.14) 



fti = fq.p = exp - 2 ^ QjXj - 2i^piyi 

3=1 1=1 



(3.15) 



On general grounds, these are eigenfunctions of the 
Laplace-Beltrami operator (and all other invariant dif- 
ferential operators) on G/K , with the eigenvalue being 



3 = 1 



1 = 1 



3=1 1=1 

-- qi{qi - 1) + 92(92 - 3) + . . . + qn{qn - 2n + I) 

Plipi + 1) - P2{P2 + 3) - . . . -Pnipn + 2n - 1), 



(3.16) 



IV. SUSY-REPLICA CORRESPONDENCE FOR 
THE n = 1 SUPERSYMMETRIC MODEL 



Let us summarize the results of two preceding sections. 
In Ref. H we reviewed Wegner's classification of poly- 
nomial spherical functions for compact replica models, 
with irreducible representations labeled by Young dia- 
grams (or sets of non-increasing positive integers giv ing 



the length of each row of the diagram). In Sec. ID 



presented an alternative classification based on the Iwa- 
sawa decomposition of the SUSY cr-model field. There, 
the A^-radial spherical functions are labeled by a set of 
non-negative integers pi and a set of parameters qj that 
are not restricted to integer or non-negative values. Ob- 
viously, the second classification is broader, in view of 
the continuous nature of the qj. Furthermore, since the 
SUSY scheme is expected to give, in some sense, a com- 
plete set of spherical functions, it should contain Weg- 
ner's classification, i.e., each Young diagram of Sec. O 
should occur as some iV-radial plane wave with a cer- 
tain set of pi and qj . We are now going to establish this 
correspondence explicitly. 

We begin with the case of minimal SUSY, n = 1. 
The starting point is a representation of Green functions 
as functional integrals over a supervector field contain- 
ing one bosonic and one fermionic component in both 
the retarded and advanced sectors. Correlation func- 
tions of bosonic fields are symmetric with respect to 
the spatial coordinates, whereas correlation functions of 
fermionic fields are antisymmetric. Thus, within the min- 
imal SUSY model one can represent correlation func- 
tions involving symmetrization over one set of variables 
and/or antisymmetrization over another set. On sim- 
ple representation-theoretic grounds, it follows that the 
n = 1 model is sufficient to make for the presence of rep- 
resentations with Young diagrams of the type shown in 
Fig. |l|. We refer to such diagrams as hooks or hook-shaped 
for obvious reasons. We introduce two "dual" notations 
(see Appendix ^ for detailed definitions) for Young di- 
agrams by counting the number of boxes either in rows 
or in columns; in the first case we put the numbers in 
round brackets and in the second case in square brack- 
ets. In particular, the hook diagram of Fig. |] is denoted 
cither as (q, 1^) or as [p+1, l'~^]. Below, we point out an 
explicit correspondence between the spherical functions 
of the n = 1 SUSY model and these hook diagrams, by 
computing the values of the quadratic Casimir operators 
and identifying the relevant physical observables. 



Evaluating the quadratic Casimir (A6) for the hook 
diagrams ((?, 1^) = [p + 1,1'~^], and taking the replica 
limit m = 0, we get 

a2((g,lP);0) =g(g-l)-2-4-...-2p 

= q{q-l)-p{p+l), (4.1) 
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-1 1 1 1 1 1 



FIG. 1. Hook-shaped Young diagrams A — (q, 1^) label the 
scaling operators that can be described within the minimal 
(n — 1) SUSY model. The numbers of boxes in each row 
and in each column are indicated to the left and above the 
diagram. The example shown in the figure corresponds to 
Q = 6, p = 4. 



which is the same (up to a constant facto r) a s the eigen- 
value of the Laplace-Beltrami operator (3^) associated 
with the plane wave ipp^q ( p.7[ ) in the SUSY formalism. 
This fully agrees with our expectations and indicates the 
required correspondence: the Young diagram of the type 
(5, 1^) of the replica formalism corresponds to the plane 
wave (pq^p (or, more precisely, to the corresponding rep- 
resentation) of the SUSY formalism. While a full proof 
of the correspondence follo ws fro m our arguments in Sec. 
VI (see especially the Sec. VIP ), we feel that the agree- 
ment between the quadratic Casimir eigenvalues is al- 
ready convincing enough for our immediate purposes. 

At this point, it is worth commenting on an apparent 
"asymmetry" between p and q in the above correspon- 
dence: the spherical function ipq^p corresponds to the 
Young diagram that has q boxes in its first row but p+1 
boxes in the first column. The reason f or th is asymmetry 
is the specific choice of positive roots (3.2). If instead of 
X — iy we chose —x + iy to be a positive root (keeping 
the other three roots), the half-sum p would change to 



p = X — ly. 



(4.2) 



This corresponds to a different choice of nilpotent sub- 
group N (generated by the root vectors corresponding 
to the positive roots) in the Iwasawa decomposition, and 
thus, to another choice of (iV-)radial coordinates x and 
y on the supcrspacc G/ K . As a result, the plane wave 



-2qx — 2ipy 



(4.3) 



characterized by quantum numbers p, q in these new co- 
ordinates, is an eigenfunction of the Laplace-Beltrami op- 
erator with eigenvalue 

i2q + if - i2p - If = Aq{q + 1) - Ap{p - 1). (4.4) 



This is the same eigenvalue as ( p.8[ ) if one makes the 
identifications q = q — 1 and p = p+1. We thus see 
that in the new coordinates the asymmetry between p 
and q is reversed: the function e~^'^~^'^^ corresponds to 
a hook Young diagram with q + 1 boxes in the first row 



and p boxes in the first column. Of course the functions 
^-2qx-2ipy g^^^^ ^-2qx-2ipy ^-^-j^ ~ ^ and p = p+1 are 

not identical (since an A'^-radial function is not iV-radial 
in general), but they belong to the same representation. 

Choosing a system of positive roots for the Iwasawa 
decomposition is just a matter of convenience; it is simply 
a choice of coordinate frame. The positive root system 
( |3.2[ ) is particularly convenient, since with this choice the 
plane waves corresponding to the most relevant operators 
(the LDOS moments) depend on x only (and not on y). 
Of course, our final results do not depend on this choice. 

We are now going to identify the physical observables 
that correspond to the operators of the n — 1 SUSY 
model. For p = the Young diagrams of the type {q, P) 
reduce to a single row with q boxes, i.e. (q) = [f], which 
in our SUSY approach represents the spherical function 
g-2qx^ As we have already mentioned, this function cor- 
responds to the moment (i^') of the LDOS. Note that 
for symmetry class A, where the global density of states 
is non-critical, the moment {v'^) has the same scaling as 
the expectation value of the q-th power of a critical wave- 
function intensity, 



Ai(r) = |^(r)| = 



(4.5) 



For the unconventional symmetry classes there is a simi- 
larly simple relation; one just has to take care of the ex- 
ponent Xp controlli ng t he scal ing of the average density 
of states, see Eqs. (LI) and (2_^). The meani ng o f the 
subscript in the notation Ai introduced in Eq. (4.5) will 
become clear momentarily. We express the equivalence 
in the scaling behavior by 



(Alir)). 
,110 



(4.6) 



We now sketch the derivatioiitll that links (i^*^) with the 
spherical function e^^'^ of the SUSY cr-model. 

The calculation of an observable (i.e. some correlation 
function of the LDOS or of wave functions) in the SUSY 
approach begins with the relevant combination of Green 
functions b eing , expressed as an integral over a super- 
vector ficld.ll3&c3 In particular, retarded and advanced 
Green functions 



Gfl,A(r,r') = iE±iTj-H)-\r,r') 



(4.7) 



(where rj is the level broadening, which for our purposes 
can be chosen to be of the order of several mean level 
spacings) are represented as 

GR{r,r')^-t{SR{r)S*j,{r')), 
GAir,r')^i{SA{r)S\ir')). (4.8) 

Here Sr^a are the bosonic components of the supcrvcctor 
field <& = {Sr, ^r, Sa, £.a) (with subscripts R, A referring 
to the retarded and advanced subspaces, respectively), 
and (. . .) on the r.h.s. of Eq. ( |4.8| ) denotes the integra- 
tion over $ with the corresponding Gaussian action of $. 
Alternatively, the Green functions can be represented by 
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using the fcrmionic (anticommuting) components ^r.a', 
wc will return to this possibility below. In order to ob- 
tain the q-th power j/' of the density of states 

'^(ro) = 7^ {Ga{vo, ro) - Gfl(ro, ro)) , (4.9) 

one has to take the corresponding combination of the 
bosonic components Si as a pre-exponential in the $ in- 
tegral: 



1 



(27r)9g: 



■( {Sniro) - 



-5^(ro))'), (4.10) 



where e*" is any unitary number. The next steps are to 
take the average over the disorder and reduce the theory 
to the non-linear cr-model form. The contractions on the 



r.h.s. of Eq. ( 4.10 ) then generate the corresponding pre- 
exponential expression in the cr-model integral: 

{i^") = 2-" {{Que - Qaa + e-'^QnA - ^'"^Qar^J , 

(4.11) 

where Q = Q(ro). The indices 6,/ refer to the boson- 
fermion decomposition. 



Although the following goes through for any value of 
a, we now take e*" ~ 1 for brevity. It is then con- 
venient to switch to Q = QA = Q(j^\ here we intro- 
duce Pauli matrices CTj in the RA space, with (T3 = A. 
It is also convenient to perform a unitary transforma- 
tion Q — ^ Q = [/ QU~^ in the RA space by the matrix 
U = {1 + iai + i<72 + ic3)/2; which cyclically permutes 
the Pauli matrices: UojU^^ = ^j-i- The combination 
of Qij entering Eq. (4.11) then becomes 



{1/2){Qrr - Qaa + Qua - QAR)bb = Qaam- (4.12) 

The Iwasawa decomposition g = nak leads to Q ~ 
na^asn'^as, where we used ka^k^^ — and aa^a^^ — 
a^cF^. Upon making the transformation Q — > Q. this 
takes the form Q 



na?a2n ^(72, or explicitly 



(l 
















\ 















1 



















1 














1 


* 

















1 













1/ 


vo 










V* 






V 



(4.13) 

where the symbol * denotes some non-zero matrix ele- 
ments of nilpotent matrices, and we have reversed the 
boson-fermion order in the advanced sector in order to 
reveal the meaning of the Iwasawa decomposition in 
the best possible way. As explained above, the vari- 
ables X and y parametrize the abelian group A which 
is non-compact in the a;-direction and compact in the y- 
direction. By observing that the 44-element of the prod- 
uct of matrices on the r.h.s. of (4.13) is e~^^, it follows 



that the matrix element ( |4.12 ) is equal to 



Qaam — e 



'2x 



(4.14) 



This completes our review of the correspondence between 
LDOS or wave- function momeats and the spherical func- 
tions of the SUSY formalism:^ 



{A\) 



V<J,0 



-2qx 



(4.15) 



Let us emphasize that, although our derivation assumes 
g to be a non- negative integer, the correspondence ( 4.15 ) 
actually holds for any complex value of q. Indeed, both 
sides of ( 4.15| ) are defined for all g £ C, and by Carl- 



son's Theorem the complex-analytic function q 1—^ iy'^) is 
uniquely determined by its values for g G N. 

At this point the unknowing reader might worry that 
the positivity of {u'') > could be in contradiction with 
the pure-scaling nature of the operator ipq^ = e~'^''^ . In- 
deed, one might argue that if the symmetry group G is 
compact, then every observable A that transforms ac- 
cording to a non-trivial irreducible representation of G 
must have zero expectation value with respect to any G- 
invariant distribution. This apparent paradox is resolved 
by observing that our symmetry group G is not compact 
(or, if fcrmionic replicas are used, that the replica trick 
is very tricky). In fact, the SUSY cr-model has a non- 
compact sector which requires regularization b y th e sec- 
ond term (or similar) in the action functional (2J_). Re- 
moving the G-symmetry breaking regularization (h — > 0) 
to evaluate observables such as (j^'), one is faced with a 
limit of the type x c» which does lead to a non-zero 
expectation value {fq,o) 7^ 0. 

We now turn to Young diagrams (P) = [p] (where we 
use the notation p = p + 1 as before) , which encode to- 
tal antisymmetrization by the permutation group. These 
correspond to the maximally antisymmetrized correla- 
tion function of wave functions. In fact, such a diagram 
gives the scaling of the expectation value of the modulus 
squared of the Slater determinant, 

Apir,,...,rp)^\Dpir,,...,rp)\^, (4.16) 
/V'i(ri) ••• Vi(rp)^ 



Dp(ri, . . . ,rp) = Det 



(4.17) 



\V'p(ri) ••• Vp(rp)y 



Here all points arc assumed to be close to each other 
(on a distance scale given by the mean free path I), so 
that after the mapping to the cr-modcl they become a sin- 
gle point. Actually, the scaling of the average (Ap) with 
system size L does not depend on the distances |r,i — rj| as 
long as all of them are kept fixed when the limit L — > 00 
is taken. However, we prefer to keep the distances suffi- 
ciently small, so that our observables reduce to local oper- 
ators of the cr-model. Moreover, all of the wave functions 
ipi are supposed to be close to each other in energy (say, 
within several level spacings). Again, larger energy dif- 
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ferences will not affect the scaling exponent; they only set 
an infrared cutoff that determines the size of the largest 
system displaying critical behavior. Clearly, Ap reduces 
to Ai = |V'(r)P when p = 1 (which was the reason for 
introducing the notation Ai above). 

In the SUSY formalism the average (Ap) can be rep- 
resented in the following way. We start with 

Ap ^ aairi) - CAiriMRiri) ~ Uri)] 

xK«(r2)-a(r2)][CK(r2)-a(r2)]... 
xKfl(r^) -a(rp)]Kfl(rp) -CA(r^)]), (4.18) 

where ^, ^* are the fermionic components of the supervec- 
tor $ used to represent electron Green functions. This 
expression can now be disorder averaged and reduced to a 
cr-model correlation function. By a calculation similar to 
that for the moment {ly'^) we now end up with the average 
of the p**^ moment of the fermion-fermion matrix element 
Qaajj of the Q matrix. Here the alternative choice of 
positive root system mentioned above (for which the ra- 
dial coordinates were denoted by a;, y) is more convenient 
since, in a sense, it interchanges the roles of x and y in 
the process of fixing the system of positive roots. As a 
result, we get the correspondence 



A, 



'Pa,p 



-2ipy 



(4.19) 



Combining the two examples above [a single-row 
Young diagram (g) and a single-column Young diagram 
(P"*"^)], one might guess that a general hook-shaped di- 
agram {q, IP) would correspond to the correlator 



{Ar^A 



p+1) 



(4.20) 



This turns out to be almost correct: the hook di agram 
(g, IP) indeed gives the leading scal ing b ehavior of (4.20). 
However, the correlation function ( 4.20| ) is in general not 
a pure scaling operator but contains subleading correc- 
tions to the scaling for {q, 1^). 

We note in this connection that in the two examples 
above each of the wave-function combinations A^ and 
Ap corresponds to a single exponential function on the cr- 
model target space, and thus to a single G-representation. 
Therefore, at the level of the a-model these combinations 
do correspond to pure scaling operators. (For the LDOS 
moments Af this was evident from the results of Ref. |l^ 
but we did not stress it there.) We will show below how 
to construct more complicated wave-function correlators 
that correspond to pure scaling operators of the cr-modcl. 



V. GENERAL WAVE-FUNCTION 
CORRELATORS 

Clearly, one can construct a variety of wave-function 
correlators that are different from the totally symmetric 
{AD and totally antisymmetric (Ap) correlators consid- 
ered in Sec. [V. One example is provided by correlation 



functions that arise when one studies the influence of ip. 
teractions on Anderson and quantum Hall transitions .E3 
In that context, one is led to consider moments of the 



Hartree-Fock matrix element (1.6), which involves the 
antisymmetrized combination (1.5) of two critical wave 
functions. In terms of the quantities Ap introduced 
above, Ref. calculated 

(A2(V'i, V'2; ri, r2)^2(V'3, i/'4; ra, r4)), (5.1) 

where the expanded notation indicates the wave func- 
tions and corresponding coordinates on which the A2 are 
constructed. Thus all four points and all four wave func- 
tions were taken to be different (although all points and 
all energies were stil l clo se to each other). To leading 
order, the correlator ( |5.lD scales in the same way as {A2) 
(where we take ipi = 'tps, ip2 = tp4,ri = ra, r2 = r4). The 
importance of the phrase "to leading order" will become 
clear in Sec. VI. 



As we discussed in Sec. IV, the scaling of the average 
{A2) is given by the representation with Young diagram 
(P) = [2]. The analysisEJ of the second moment {A2) 
shows that its leading behavior is given by the diagram 
(2^) = [2^]. A natural generalization of this is the follow- 
ing proposition: the Young diagram 



A [Pl,P2, ■■■,Pn 



(5.2) 



relates to the replica cr-model operator that describes the 
leading scaling behavior of the correlation function 



{Ap-^ Ap^ ■ ■ ■ Ap^ 



(5.3) 



We will argue in Sec. VI below that this is indeed correct. 
Here we wish to add a few comments. 

In general, all combinations Ap. may contain different 
points and different wave functions (as long as the points 
and the energies are close) without changing the leading 
scaling behavior. Thus, a general correlator correspond- 
ing to a Young diagram A will involve |A| points and the 
same number of wave functions. However, if 



A = bi,P2, 



(5.4) 



we may choose to use the same points and wave functions 
for all Qi combinations Ak^ of a given size ki . This yields 
a somewhat simpler correlator 



(5.5) 



with the same leading scaling. If we use the alternative 
notation 



X = iquq2,...,qn) = il1\...,l''/) 



(5.6) 



for the Young diagram (5.2), the correlator (5.5) can also 
be written as 



-1-62 ' ■ '^fci-l-... -1-6,-1^61 -I-...-I-6J' \^-') 



11 



see Eqs. (A2) and (A5) in Appendix^. In fact, as is easy 
to see, this can also be rewritten in a natural way as 



If we introduce the notation 

A, 



vi = Ai, 



A., 



2<i<n, 



(5.8) 



(5.9) 



then the correlator K\ can also be cast in the following 
form: 



K, 



(5.10) 



Below we will establish the correspondence of the cor- 
relation functions (5.10) with Young diagrams that was 
stated in this section. We will also show how to build 
pure-scaling correlation functions and establish a con- 
nection with the Fourier analysis on the symmetric space 
of the SUSY cr-model. 



VI. EXACT SCALING OPERATORS 

Let us now come back to the issue of exact scaling 
operato rs. In the preceding section we wrote down a large 
family (5.3) of wave-function correlators. In general, the 
members of this family do not show pure scaling. We are 
now going to argue, however, that if we appropriately 
symmetrize (or appropriately choose) the points or wave 
functions that enter the correlation function, then pure 
power-law scaling does hold. 



A. An example 

Let us begin with the simplest example illustrating the 
fact stated above. This example is worked out in detail in 
Sec. 3.3.3 of Ref. Q and is provided by the correlation 
function 



(|^i(ri)V2(r2)n. 



.1) 



When the two points and the two wave functions are 
different, this yields 

{{QhRM — QaAm)'^) = 2 — '2{Q RRmQ AA.bb) (6-2) 

after the transformation to the cr-model. Now the expres- 
sion 1 — Q RRfibQ AA,bb is not a pure-scaling a-model op- 
erator: by decomposing it according to representations, 
one finds that it contains not only the leading term with 
Young diagram (2), but also the subleading one, (1,1). 
To get the exact scaling operator for (2), which is 

1 — {Q RR,bbQ AA,bb + QrAmQaRm), (6.3) 

one has to symmetrize the product of wave functions in 



the correlator that does exhibit pure scaling is 

(|^i(ri)V'2(r2) + M^2)M^i)n- (6.4) 

Alternatively, one can take the points to be equal and 
consider the correlation function 



{\Mri)Mri)\'). 



(6.5) 



Then one gets the exact scaling operator right away, since 
the correlation function ( |6.5[ ) already has the required 
symmetry. One can also take the same wave function: 

(|^i(ri)^i(r2)P), (iVi(ri)O. (6.6) 

All of these reduce to the same exact scaling operator 
.31) in the cr-modcl approximation. 



B. Statement of result 

The example above gives us a good indication of how 
to get wave-function correlators corresponding to pure- 
scaling operators: the product of wave functions should 
be appropriately (anti) symmetrized before the square of 
the absolute value is taken. To be precise, in order to get 
a pure-scaling correlation function for the diagram ( |5.2| ) 
[giving the leading scaling contribution to Eq. (5.3)], one 
should proceed in the following way: 

(i) View the points and wave functions as filling the 
Young diagram (5^) by forming the normal Young 
tableau Tq (see Appendix ^ for definitions); 

(ii) Consider the product of wave-function amplitudes 



'4'i{ri)ip2(r2) ■ ■ -ipNirN), N ^ pi + . . . + p,n. (6.7) 

In the notation of Appendix |^ this is 'i' x{TQ,Tf)). 

(iii) Perform the Young symmetrization c\ = h\a\ ac- 
cording to the rules described in Appendix ^ (sym- 
metrization ax with respect to all points in each row 
followed by antisymmetrization 6a with respect to 
all points in each column). In this way we obtain 



*A(To,CAro). 



(6.8) 



(iv) Take the absolute value squared of the resulting 
expression: 



|*a(To,caTo)| 



(6.9) 



(S.l) with respect to points (or wave function indices) 



Several comments are in order here. First, one can 
define several slightly different procedures of Young sym- 
metrization. Specifically, one can perform it with respect 
to points (as described above) or, alternatively, with 
respect to wave functions (obtaining |vI'a(caTo, To)p). 
Also, one can perform the Young symmetrization in 
the opposite order (first antisymmetrization along the 
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columns, then symmetrization along the rows: c\ = 
a\b\). In fact, it is not difficult to see that carrying 
out c\ with respect to wave functions is the same as 
performin g c\ w ith r espect to points, and vice versa, 
see Eqs. (A24), (A25). While for different schemes one 



will in general obtain from (6^) slightly different expres- 
sions, they will scale in the same way upon averaging, 
as they belong to the same irreducible representation. 
Furth erm ore, once a Young symmetrization of the pro! 
duct (|6.7|) has been performed, one can, instead of tak- 
ing the absolute value squared, simply multiply it with 
the product ?/;5'(ri)V'2(r2) • ■ • V'Ar(rjv)- Finally, the sym- 
metrization with respect to points is redundant if the 
corresponding poin ts (o r wave functions) are taken to be 
the same, see Eq. (A2S). 

To illustrate the procedure, let us return again to the 
correlation function (^.1[) 



(6.10) 



(^2 ( V'l , V'2 ; r 1 , r2 ) ^2 (V'a , -04 ; rs , r4 ) ) , 



considered in Ref. [Q. As we have already discussed, its 
leading scaling is that of the Young diagram (2^); how- 
ever, Eq. ( 6.10| ) includes also corrections due to sublead- 
ing operators. In order to get the corresponding pure- 
scaling correlation function, we should start from the 
product ^i(ri)?A2(r2)'/'3(r3)V'4(r4) and apply the Young 
symmetrization rules corresponding to the diagram (2^). 
This will lead to the expression 



[V'l(l"l)V'2(r2) + '0l(l"2)'02(ri)] 

X [V'3(r3)V'4(r4) + V3(r4)V'4(r3)], 



.11) 



further anti-symmetrized with respect to the interchange 
of ri with ra and with respect to interchange of r2 with 
r4. Finally, one should take the absolute value squared. 
As an alternative to the symmetrization, one can simply 
set ri = r2 and = r4 (which means choosing the min- 
imal Young tableau Tmin for T^), in which case there is 
no need to symmetrize. This results in 



|[^i(ri)V'3(r3) - ^i(r3)^3(ri)] 
X [■02(ri)'04(r3) - V'2(r3)V'4(ri)]| 



3.12) 



A similar expression can be gotten by setting ipi = 'ip2 
and ■i/'3 = ^4- Finally, one can do both, keeping only two 
points and two wave functions. This results exactly in 



^'(22) (Tmin, b(2^)T^ 



= \DiV ^A. 



2i 



3.13) 



which is thus a pure scaling operator. 

This has a natural gen eralization to the higher-order 
correlation functions (5.8) as follows. Let ri, . . . , r„ be a 
set of n distinct points. For each m < n evaluate A,„ at 
the point r„i on a set of wave functions • ■ • j V'rn"''- 

The coincidence of evaluation points takes care of the 
symmetrization along all rows of the Young diagram. 
Moreover, the antisymmetrization is included in the def- 
inition of the Ai. Therefore, with such a choice of points 



the correlation function (5.8) will show pure scaling. This 
statement is independent of the choice of wave functions 
V'j'"'': all of them can be different, or some of them cor- 
responding to different m can be taken to be equal. The 
most "economical" choice is to take only n different wave 



functions tpi, . . . ^tpn and for each m set ip 
pendent of m) for j = 1, . 



(m) 



i/jj (inde- 



3 -rj 

m. This is the choice made 



by the minimal tableau (see Eq. (A29)): 



^91^92 ^92 -93 



3.14) 



where the numbers (gi, . . . , q„) specify the representation 
with Young diagram A as in Eq. (|5.6|). 



C. Sketch of proof 

We now sketch the proof of the relation between the 
wave-function correlation functions with the proper sym- 
metry and the cr-model operators from th e cor responding 
representation. In accordance with Eq. (4.10), we begin 
with an integral over a supervcctor field S, 



{cx{S^{r,)---S]^irN)}cx{Sr{ri 



5* 



N (rjv)}). 

(6.15) 

As before, S denotes the bosonic components of the su- 
perfield; the superscript in S~ reflects the structure in the 
advanced-retarded space: S~ = Sr — Sa- This structure 
ensures that, upon performing contractions, we get the 
required combinations of Green functions, Gr — Ga- We 
emphasize, however, that one could equally well choose 
5g + Sa or 5*^ — e*"5A for any a, as was done in Eq. 
( |4.10| ). Indeed, by Eq. (|||) all that matters is that the 



coefficients of Sr and Sa have the same absolute value. 
We also mention that the freedom in choosing a is eluci- 
dated in more detail in Sec. VI F and Appendix B 2. 

The subscript of the S'-fields in Eq. ( |6.15D is the replica 
index. (Recall that we consider an enlarged number 
of field components.) The symbol cx{. . .} denotes the 
Young symmetrization of the replica indices according to 
the chosen Young diagram A = (gi, 52, • ■ ■) = [pitP2t ■ ], 
and N = \\\ = Yl, Pi — Ij ■ It is given by the prod- 
uct cx = bxax of the corresponding symmetrizatio n and 
antisymmetrization operators. (Although in Eq. (6.15) 
we put Cx twice, it would actually be sufficient to Young- 
symmetrize only S fields, or only S* fie lds.) It is possi- 



ble to express the correlation function ( 3.15 ) in a more 
economical way (i.e., by introducing fewer field compo- 
nents), without changing the scaling operator that results 
on passing to the cr-model. This economy of description 
is achieved by observing that symmetrization is provided 
simply by the repeated use of the same replica index: 

{bx{S^{r['^) ■ ■ ■ 5r(rW)52-(rf )) • • • S^ (r^^y) ■ ■ ■ 



xS-{r["^)---S-{ri:'>)} 

xbx{Sr{r['^) ■ ■ ■ (rW)52*-(rf )) • • • S^r^^J) 

x5r(ri"^)---5r(r[::^)})- (e.ie) 
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Here we denoted by r'-''' , . . . , Tq-'^^ the points filling the j- 
th row of the Young diagram (gi, . . . g„) = [pi, • ■ • ,Pm], 
and b\{. . .} still denotes the operation of antisymmetriza- 
tion along the columns of the Young diagram. Perform- 
ing all Wick contractions and writing Green fu nctio ns as 
s ums over wave functions, one sees that Eqs. (6.15) and 
( |6.16 ) give (up to an irrelevant overall factor) exactly the 
Young-symmetrized correlation function of wave func- 
tions that was described in Sec. 0. Specifically, the ob- 
tained correlation function yields the average of Eq. (|6.9[). 

By the process of transfo rming t o the cr-model, the 27V 
field values of S and S* in ( 6.15 ), ( 6. IE ) get paired up in 
all possible ways to form a polynomial of 7V-th orderj-in 
the matrix elements of Q. The general rule for this isE3 



is a pure-scaling operator of the non-linear cr-model. This 
will be achieved by showing that Pa is an eigenfunction 
of all Laplace-Casimir operators for G/K. The latter 
can be done in two different ways. Firstly, one may ar- 
gue with the help of the Iwasawa decomposition that Px 
is an A^-radial spherical function and thus has the de- 
sired eigenfunction property. In subsection VXD below. 



we spell out this argument along with its natural gener- 
alization to complex powers q. Secondly, it is possible to 
get the desired result directly (without invoking the Iwa- 
sawa decomposition) by showing that the function P\ is a 
highest-weight vector for the action of G on the matrices 



Q. This is done in subsection VI F 



S-ir,)S;;ir2) ^ /(|ri -r2|)Qp,p, {'2(^1+^2)) , 

(6.17) 

where the prefactor /(|ri — r2|) = {ttv) ^Im(Gyi(ri, r2)) 
depends on the distance between the two points, and 

Q = Qaam = hi QBB. - Qaa + Qra - QAR)bb was in- 
troduced in Eq. ( 4.12 ). In a 2D system for example, 
/(r) = e~^^^''Jo{kpr). The key properties of the function 
/(r) are /(O) = 1 (more generally, /(r) ~ 1 as long as 
the distance is much smaller than the Fermi wave length, 
r <^ Xf) and /(r) <C 1 for r ^ A_f. In the latter case the 
corresponding pairing between the fields S and S* can 
be neglected. Assuming that all points in the correlation 
function (5.15) are separated by distances r 3> Ai?, we 



get an expression of the diagonal structure 



){QiiQ2 



iNN 



(6.18) 



where c 



means that we Young-symmetrize sepa- 



rately wit h resp ect to both sets of indices (left and right). 
(If in Eq. ( |6T5| ) only one Young symmetrizer is included, 
then only the corresponding set of indices is Young sym- 
metrized here; this does not change the irreducible repre- 
sentation that Eq. (6TS) belongs to.) Similarly, starting 
from Eq. (3.16) and assuming that all points are suffi- 
ciently well separated, we obtain 



((4^'®4«^)(Q,,,,Q,,,,...Q,„,„)), 



3.19) 



where the first qi indices ji are equal to 1, the next q2 
are equal to 2, and so on, and the last q„ are equal to n. 
In this case the operator a\ for symmetrization is redun- 
dant (since symmetrization of equal indices has a trivial 
effect) and we may simplify the expression by replacing 
c\ by the operator 6a for antisymmetrization along the 
columns of the Young diagram. O ne ca n also t ake some 
points in the original expressions ( 6.15| ), ( |6.16 ) to coin- 
cide (provided that the result does not vanish upon an- 
tisymmetrization); this will not influence the symmetry 
and scaling nature of the resulting correlation functions. 

To complete our (sketch of) proof, we must show that 
the polynomial 



Pa 



(L) 
A 



Qi.j.) (6.20) 



D. Argument via Iwasawa decomposition 

We now argue that the polynomial Pa is an eigenfunc- 
tion of all Laplace-Casimir operators for G/K. To this 
end, our key observation is that Pa can be written as a 
product of powers of the principal minors (i.e., in our 
case, the determinants of the right lower square sub- 
matrices) of the matrix Q = QAA,bb for the case of n 



replic as. Indeed, following the derivation of Eqs. ( 3.14 ), 
( A2£ ) , we can associate the left indices of the n x n ma- 
trix Q with one minimal Young tableau, and the right 
indices with another minimal tableau. As a result, if we 
denote by dj the principal minor of Q of size j x j, we 
see that 



q2 



df-^'-^dir, 



(6.21) 



since the Young symmetrizer ca here acts essentially as 
the antisymmctrizcr 6a, producing determinants of the 
principal submatrices of Q. 

The final step of the argument is to show that Pa agrees 
(up to a co nstan t) with the A^-radial spherical function 
^q^o of Eq. §J^: 



(6.22) 



which is already known to have the desired property. For 
that, let us write Q in Iwasawa decomposition as 



/l 





1 * 
ly 



-2x„ 







\ / 





V 



-,-22:2 



e — 




7 



o\ 



1 



= 1/ 

(6.23) 



(Precisely speaking, this is the Iwasawa decomposition of 
the full matrix Q = ha? a2fi^^ <J2 projected t o the boson- 
boson part of the right lower block; cf. Eq. ( 4.13| ).) Due 
to the triangular form of the first and last matrices in this 
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decomposition, the principal minors dj of this matrix are 



-2xi 



= exp 



3.24) 



When this expression is substituted into Eq. (6.21), we 
get exactly the function tp^ ^ of Eq. (3.15) for the set q = 
(gi, . . . , g„) of positive integers qj. Since this function 
is an eigenfunction of all Laplace-Casimir operators for 
G/K, it follows that P\ has the same property. This 
completes our proof. 



To summarize, recall that in Sec. VI E we specified a 
certain set of wave-function correlators. Our achievement 
here is that we have related these correlators to pure- 
scaling operators of the non-linear cr-model. By doing so, 
we have arrived at the prediction that our wave-function 
correlators exhibit the same pure-power scaling. 

Finally, let us remark that, although the analysis above 
was formulated in the language of the SUSY cr-model, 
it could have been done equally well for the replica a- 
models. (In the presence of a compact sector, where 
the Iwasawa decomposition is not available without com- 
plexification, it would actually be more appropriate to 
carry out the final step of the argument by th e theory 
of highest- weight vector as outlined in Sec. VI F and Ap- 
pendix pi) 



definition of and the fact that the total density of states 
is self-averaging, one has an a-priori bound for positive 
real values of q : 



< f{q) < {L/lfi [q > 0), 



(6.27) 



where I is the lattice spacing (or UV cutoff) of the d- 
dimensijOiial system. In conjunction with the functional 
relatiorO f{q) = f{l—q), this inequality leads to a bound 
of the form 



|/(g)|<e^-(i+|i^''«l) (qeC), 



(6.28) 



where Cl oc InL is a constant. Thus, in finite volume, 
/ is an entire function of exponential type and is also 
bounded along the imaginary axis. By Carlson's The- 
orem, this implies that / is uniquely determined by its 
values on the non-negative integers. It follows that the 
result of our derivation, taking the pure-scaling correla- 
tion functions ( ^.8D to cr-model expe ctatio n values of the 
A'^-radial spherical functions ipq^ of (3.25), extends from 
non-negative integer values of q to all complex values of q. 
This relation is expected to persist in the infinite- volume 
limit i — > oo. 



Alternative construction of scaling operators: 
highest-weight vectors 



E. Generalization to arbitrary qj 

We now come to a generalization of our correspon- 
dence. The correlators considered in Sec. VI up to now 
were polynomials (of even order) in wave-function am- 
plitudes ip and ijj* , and the resulting cr-model operators 
were polynomials in Q. The important point to empha- 
size here is that the wave-function correlation functions 



(5.8) are perfectly well-defined for all complex values of 
the exponents qj (j ~ I, . . . ,n). At the same time, while 
the polynomial a-model operators of Wegner's classifica- 
tion clearly require the numbers qj to be non-negative 
integers, the A''-radial spherical functions (3.15) given by 
the SUSY formalism. 



ipqfi = exp 



.1=1 



3.25) 



do exist for arbitrary quantum numbers q = (gi, . . . , g„). 
Thus one may suspect that our correspondence extends 
beyond the integers to all values of q. This turns out to 
be true by uniqueness of analytic continuation, as follows. 



We gave an indication of the argument in Sec. [V and 
will now provide more detail. Let n = 1 for simplicity 
(the reasoning for higher n is no different), and consider 

f{q) EE {u'^)/{vY {q e C). (6.26) 

The triangle inequality gives \f{q)\ < /(Reg). By the 



In previous sections we constructed scaling operators 
in the cr-model from the Iwasawa decomposition. Here 
we show how to construct the same operators by using a 
different approach based on the notion of highest-weight 
vector. We just outline the basic idea of this approach, 
relegating details of the construction to Appendix 

The cr-model field Q takes values in a symmetric space 
G/K. Our goal is to identify gradientless scaling oper- 
ators of the CT- model, i.e., operators that reproduce (up 
to multiplication by a constant) under transformations 
of the renormalization group. We know that the change 
of a local cr-model operator, say A, under an infinitesi- 
mal RG transformation can be expressed by differential 
operators acting on A considered as a function on G/K. 
Assuming that the cr-model Lagrangian is G-invariant, 
the infinitesimal RG action is by differential operators 
which are G-invariant (also known as Laplace-Casimir 
operators). Thus a gradientless operator of the cr-model 
is a pure scaling operator if it is an eigenfunction of the 
full set of Laplace-Casimir operators on G/K . 

Such cigenfunctions can be constructed by exploiting 
the notion of highest- weight vector, as follows. Let g = 
Qc denote the complexified Lie algebra of the Lie group 
G. The elements X € Q act on functions f{Q) on G/K 
as first-order differential operators X: 



{Xf){Q) = 



dt 



(6.29) 



By definition, this action preserves the commutation re- 



15 



lations: [X,Y] ^ [X,Y]. 

Fixing a Cartan subalgebra t) C g we get a root-space 
decomposition 



= n+ e f) I 



(6.30) 



where the nilpotent Lie algebras n± are generated by 
positive and negative root vectors. We refer to elements 
of n+ (n_) as raising (resp. lowering) operators. (Com- 



paring with the Iwasawa decomposition of Sec. Ill, we 
observe that n+ is the same as the complcxification of n, 
and a is a subalgebra of f), with the additional generators 
of f) lying in the complexified Lie algebra of K .) 

Now suppose that ip\ is a function on G/K with the 
properties 

1. Xlpx for all X e n+, 

2. Hlpx = \{H)lpx for ah e t (6.31) 

Thus (fx is annihilated by the raising operators from n+ 
and is an eigenfunction of the Cartan generators from fj. 
Such an object ip\ is called a highest-weight vector, and 
the eigenvalue A is called a highest weight. 

Since the Lie algebra acts on functions on G/K by 
first-order differential operators, it immediately follows 
that the product ip\^+\^ = ipx^tpx^ of two highest- weight 
vectors, as well as an arbitrary power (pqx = of a 
highest- weight vector, are again highest- weight vectors 
with highest weights Ai + A2 and gA, respectively. In the 
compact case the power q has to be quantized (a non- 
negative integer) so that ip\ is defined globally on the 
space G/K. On the other hand, in the non-compact case, 
we can find a positive ((/?a > 0) highest-weight vector, and 
then it can be raised to an arbitrary complex power q. 

Now recall that a Casimir invariant C is a polynomial 
in the generators of g with the property that [C, X] = 
for all X e g. The Laplace-Casimir operator C is the 
invariant differential operator which corresponds to the 



Casimir invariant C by the action ( 6.29 ). If a function ^px 
has the highest- weight properties ( p. 31 ), then this func- 
tion is an eigenfunction of all Laplacc-Casimir operators 
of G. To sec this, one observes that on general grounds 
every Casimir invariant G can be expressed as 



G — G^f + DaXa, 



(6.32) 



a>0 



where every summand in the second term on the right- 
hand side contains sonic Xa G n+ as a right factor. Thus 
the second term annihilates the highest- weight vector ipx ■ 
The first term, C7|,, is a polynomial in the generators of 
the commutative algebra f) and thu s has ipx as an eigen- 
function by the second relation in (6.31). 

In summary, gradientless scaling operators can be 
constructed as functions that have the properties of a 
highest-weight vector. To generate the whole set of such 
operators, one uses the fact that powers and products 
of hcighcst-weight vectors are again hcighest-wcight vec- 



tors. 

Let us discuss how this construction is related to the 
Iwasawa decomposition G = NAK. A^-radial functions 
f{Q) on G/K by definition have the invariancc property 



/(nQn-i) = f{Q) 



\fn e TV. 



(6.33) 



Any such function is automatically a highest-weight vec- 
tor if the nilpotent group N is such that its (complexified) 
Lie algebra coincides with the algebra n+ of raising op- 
erators. Indeed, if X is an element of the Lie algebra of 
N, then 



iXf)iQ) = - 



0, 



(6.34) 



since the expression under t he t d erivative does not de- 
pend on t by the invariance ( 6.33| ). 

In Appendix |b| we implement this construction explic- 
itly. We consider certain linear functions of the matrix 
elements of Q, which we write as 

fiY{Q)^TT{YQ). (6.35) 

From the definition ( 6.29| ) it is easy to see that 



{XfiY){Q) = ^ ^^Jr(e*^ye-*^0) = M[x.F](Q)- 

(6.36) 

Then, if [v'T, F] = 0, the function /iy(Q) is annihilated 
by X. To construct highest-weight vectors, which are 
annihilated by all X for X G n+, we then build certain 
polynomials from these linear functions, and form prod- 
ucts of their powers. In this man n er we recover exactly 
the set of scahng operators ( 3.21 ), ( 6.25 ). 



VII. WEYL GROUP AND SYMMETRY 
RELATIONS BETWEEN SCALING EXPONENTS 

In the preceding sections we constructed wave- function 
observables that show pure-power scaling, by establish- 
ing their correspondence with scaling operators of the 
SUSY (T-model. Now we are ready to explore the im- 
pact of Weyl-group invariance on the spectrum of scaling 
exponents for these operators (and the corresponding ob- 
servables) at criticality. The Weyl group W is a discrete 
group acting on the Lie algebra a of the group A, or 
equivalently, on its dual a*. Acting on a*, W is gener- 
ated by reflections Va at the hyperplanes orthogonal to 
the even roots a: 



/i I— > yU — 2a 



[a, a) ' 



(7.1) 



where (•, •) is the Euclidean scalar product of the Eu- 
clidean vector space a*. 

Key to the following is the Harish-Chandra isomor- 
phism, see Refs. [pSl p6i for the classical version and 
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Ref. 1 28 for the SUSY generalization (which we need). 
The statement is that there exists a homomorphism (ac- 
tually, an isomorphism in classical situations) from the al- 
gebra of G-invariant differential operators on G/K to the 
algebra of M^-invariant differential operators on A. This 
homomorphism (or isomorphism, as the case may be) is 
easy to describe: given a G-invariant differential operator 
D on G/K, one restricts it to its iV-radial part, which 
can be viewed as a differential operator on A, and then 
performs a so-called Harish-Chandra shift (A — >■ \ — p) 
by the half-sum of positive roots p. The shifted operator 
turns out to be M^-invariant. 

This property of H^-invariance is what matters to us 
here, for it has the consequence that if XtJ.{D) denotes 
the eigenvalue of D on the s pher ical function (or highest- 
weight vector) (ys^, see Eq. (3.15), then 



this general statement. As bef ore, w e focus on class A 



(7.2) 



for all w € W. In words: if two spherical functions ipi_i 
and ipx have highest weights A = wp related by a Weyl- 
group element w G VF, then their eigenvalues are the 
same, Xm(^) ~ X\{D), for any D. To the extent that 
the cr-model renormalization group transformation is G- 
invariant (and hence is generated by some G-invariant 
differential operator on G/K), we have the following 
important consequence: the scaling dimensions of the 
scaling operators (which arise as eigenvalues of the G- 
invariant operator associated with the fixed point of the 
RG flow) are W^-invariant. 

For our purposes it will be sufficient to focus on the 
subgroup of the Weyl group which is generated by the 
following transformations on a* : (i) sign inversion of any 



one of the /i^-components: 



-p!l (reflection at the 



hyperplane p^ = 0), and (ii) pairwise exchange of pP- 
components: /i? -s-> /i^ (reflection at the hyperplane p^ — 



p, = 0). In view of Eq. (3.14) these induce the following 



transformations of the plane- wave numbers qj : 

(i) sign inversion of qj + for any j e {1, 2, . . . , n}: 



(7.3) 



where Cj is the coefficient in front of Xj in the ex- 
pressi on fo r the h alf-sum p of positive roots, see 
Eqs. (|Ill), (||l|); 



(ii) permutation of q^ 
i,j e {l,2,...,n}: 



and 



1j 



for some pair 



(7.4) 



By combining all such operations, one generates a sub- 
group Wo C of the Weyl group. Whenever two scaling 
operators with quantum numbers q = (gi, g2, ■ • ■ , Qn) and 
q' = {q[, q2, ■ . ■ , q'j-^) are related by a Weyl transformation 
w G Wq, the scaling dimensions of these scaling opera- 
tors must be equal. We now present some examples of 



for which Cj = 1 — 2j, see Eq. (3.12). Generalizations to 



the other classes will be discussed below. 

Consider first the most symmetric representations (q), 
which are characterized by a single number qi = q. (Here, 
for convenience, we continue to use Young-diagram no- 
tation, even though qi need not be a positive integer and 
does not correspond to a representation of polynomial 
type.) The invariance under the Weyl group then im- 
plies that the following two representations 



(9), (1-g), 



(7.5) 



(here we used ci = —1) give identical scali ng d imensions. 
This is exactly the symmetry statement ( |l.4| ) governing 
the multifractal scaling of the LDOS moments. 

Next, consider representations of the type (gi, 92)- By 
applying the Weyl symmetry operations above, we can 
generate from it a series of 8 representations: 

(91, 92), (1-91,172), (91,3-^2), (1-91,3-52), 

(2-52,2-91), (-1 + 92,2-91), 

(2-92,1 + 91), (-1 + 92,1 + 91)- (7.6) 

Again, all of them are predicted to give the same scaling 
dimension. As an important example, starting from the 
trivial representation (0) = (0, 0) (i.e. the unit operator) 
we generate the following set: 



(0,0) 
(2,2) 



(1,0), (0,3), 
(-1,2), (2,1) 



(1,3), 
> (-1,1). 



(7.7) 



Since (0, 0) has scaling dimension zero, we expect the 
same to hold for all other representations of this list - as 
long as the Anderson-transition fixed point is in class A. 
This is a remarkable statement. 



In fact, among the set of representations (7.7), four 



arc of polynomial type and "standard" in that they are 
also present in the replica approach of Wegner. Aside 
from the trivial representation (0,0), these are (1,0), 
(2,2), and (2,1). The representation (1,0) corresponds 
to (Q) which is well known to be non-critical in the 
replica limit. However, for the polynomial representa- 
tions (2, 2) and (2, 1), our exact result seems to be new. It 
is worth emphasizing that Wegner's four-loop perturba- 
tive C-functioncJ is fully consistent with our finding: both 
coefficients 02 and C3 vanish for these operators in the 
replica Uyit, see Table || above. Moreover, a numerical 
analysifJ of the correlation function (5.1), whose lead- 
ing scaling behavior is controlled by (2, 2), also yielded a 
result consistent with X(2.2) = 0. 

For our next example of imp ortance, consider the case 
of 9i = 92. Inspecting Eq. ( [7.6[ ) we see, in particular, that 
the cr-model operators (9, 9) and (2 — 9, 2 — 9) have the 
same sca ling dimensions. Now we know [see Sec. VI D 
and Eq. ( |5.8| )] that the operator (9, 9) corresponds to the 
moment of t he Hartree-Fock type correlation func- 
tion A2, Eq. ( |l.5| ). Thus we learn that the multifractal 
spectrum of scaling dimensions for the Hartree-Fock mo- 
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ments A2 is symmetric under the reflection go 2 — g. 

One can continue these considerations and look at 
equivalences between scaling dimensions for n = 3, i.e. 
for operators {qi, (72, ^s), and so on. In general, the Weyl 
orbit of an operator (qi, . . . , qn) with n different compo- 
nents consists of 2"ri,! operators (due to 2" sign inversions 
and n\ permutations) with equal scaling dimcnsiopa. Wc 
have checked that all results obtained by Wcgncr,Ej who 
analyzed operators described by Young diagrams up to 
size 5 and up to four loops (see Table ^ , are fully consis- 
tent with this prediction. 

This includes the afore- mentioned representations (1), 
(2, 1), and (2, 2), as well as (3, 1, 1): all of them are re- 
lated to the trivial representation by Weyl-group opera- 
tions and, indeed, Wegner obtained zero values of 02 and 
C3 for all of them (in the replica limit). Furthermore, the 
operators (3, 2) = [2, 2, 1] and (3, 1) = [2, 1, 1] are clearly 
related to each other by the Weyl reflection (72 — >■ 3 — (72 . 
Again, as is shown in Table |, Wegner's four-loop results, 
a2 — A and C3 — 24, are the same for these. 



VIII. OTHER SYMMETRY CLASSES 

In order to apply the Weyl-symmetry argument to the 
other symmetry classes, we need the expressions for the 
half-sum of positive roots for them. More specifically, we 
will now present the "bosonic" part pb (which is a linear 
combination of the basic functions Xj) of p. By transcrip- 
tion of the above analysis, its coefficients Cj determine 
the Harish-Chandra shift entering the Weyl transforma- 
tion rules for the operators (gi, . . . , g„), see Eqs. (7.5) 
and (|7j). 

Th e ro ot systems for all symmetry classes arc listed in 
Table III of Appendix The resulting pf, are 



Pb 



(8.1) 



where the coefficients Cj {j ~ 1,2, 



read 





= 1 - 2j, 


class A, 


(8.2) 




= -j, 


class AI, 


(8.3) 


Cj 


= 3 - 4j, 


class All, 


(8.4) 


Cj 


= 1 - 4j, 


class C, 


(8.5) 


Cj 


= l-j, 


class D, 


(8.6) 


Cj 


= -2j, 


class CI, 


(8.7) 


Cj 


= 2 - 2j, 


class Dili, 


(8.8) 




1 

- 2 ~ J' 


class BDI, 


(8.9) 


Cj 


- 2 - 4j, 


class CII, 


(8.10) 




= 1 - 2j, 


class AIII. 


(8.11) 



The results obtained above for class A generalize in a 
straightforward manner to four of the other classes, which 
comprise the two remaining Wigncr-Dyson classes, AI 
and All, and two of the Bogoliubov-de Gcnnes classes. 



C and CI. The Weyl-symmetry operations involve the 
pertinent values of cj in each case. For example, for the 
most symmetric operators (g) (characterizing the LDOS 
moments) we obtain the correspondence (g) (— ci — g), 
where — ci has value 1 for the classes A, AI, and All, value 
2 for class CI, a nd value 3 for class C. This is exactly the 
symmetry ( |1.4| ) obtained in Rcf. with g* — ~c\. 

Correspondences between the representations with two 
or more numbers (gi, . . . , g„) are obtained in exactly the 
same way as described in Sec. VII for class A. Agai», 
we have checked that the four-loop results of WegnerE2l 
for the orthogonal and symplectic classes (AI and All) 
conform with our exact symmetry relations. Specifi- 
cally, for class AI, our results imply the following Weyl- 
symmetry relations (and thus equal values of the scaling 
dimensions): (i) (2,2) ^ (2); (ii) (1,1) o (2,1,1); (iii) 
(3,2) o (3); (iv) (2,2,1) o (1) O (0) (scaling exponent 
equal to zero); these are the Young diagrams up to size 
5 studied in Rcf. |2^. For class All the dual correspon- 
dences hold: (i) (2,2) o (1,1); (ii) (2) o (3,1); (iii) 
(2,2,1) o (1,1,1); (iv) (3,2) o (1) ^ (0). Needless to 
say, the results of Ref. for the coefficients 02 , C3 for 
these operators do conform with the predicted relations. 

Generalization to the remaining five classes (D, Dili, 
BDI, CII, and AIII) is more subtle due to peculiarities of 
their cr-model manifolds. We defer this issue to Sec. p3. 



IX. TRANSPORT OBSERVABLES 

We now address the question whether the classification 
and symmetry analysis of the present paper are also re- 
flected in transport observables. To begin, we remind the 
reader that such a correspondence between wave- function 
and transport observables has previously been found for 
the case of the (g) operators. Specifically, one can con- 
sider the scaling oL-ppmcnts of the two-point conduc- 
tance at criticalityj^l^ 



(.9'(r,r')) 



r r 



(9.1) 



Actually, g'^{v, r') is not a pure-scaling operatoicJ (u nlike 
the LDOS moments considered above), thus Eq. (9.1) 
should be understood as characterizing the leading long- 
distance behavior of (g'^(r,r')). Nevertheless, it turned 
out that the transport exponents F^ j^nd the LDOS ex- 
ponents Xq — ^q^ qxp are related asl 




g < g*/2, 
g > g*/2. 



(9.2) 



Notice that while the LDOS spectrum Xq is symmetric 
with respect to the point g*/2 = — ci/2, the two-point 
conductance spectrum F^ "terminates" (i.e., has a non- 
analyticity and becomes constant) at this point. Yet, 
the spectrum Tq clearly carries information about the 
Weyl symmetry: if one performs its analytic continua- 
tion (starting from the region below g«/2), one gets the 
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spectrum 2xq = 2xq*-q. 

A physically intuitive argument explaining Eq. ( |9.2[ ) is 
as follows. For sufficiently low the moments of (/(r, r') 
are controlled by small values of the conductance. When 
^(r, r') is small, one can think of it as a tunneling conduc- 
tance that is proportional to the product of the LDOS 
at the points r and r'. The corresponding correlation 
function {v'^{v)v'^[y')) scales with |r — r'| with the e xpo - 



where 



nent 2xq^ in agreement with the first line of Eq. (3.2). 



(This argument ca n be also cast in the RG language, 
see the end of Sec. IX B.) On the other hand, the two- 



point conductance cannot be larger than unity. For this 
reason the relation = 2xq does not hold beyond the 
symmetry point q = The moments with q > (7*/2 

are controlled by the probability to have g{r, r') of order 
unity. 



In view of the relation (9.2), a natural question is 
whether there are any transport observables correspond- 
ing to the composite operators (gi,g2,---) beyond the 
dominant one, (g). We argue below that this is indeed 
the case, construct explicitly these transport observables, 
and conjecture a relation between the critical exponents. 

In order to get some insight into this problem, it is 
instructive to look first at quasi-lD metallic systems, 
whose transport pimaerties can be described within the 
DMPK formalism.ElO The rationale behind this is as fol- 
lows. First, the classification of transport observables 
that we are aiming at is based (in analogy with the 
classification of wave-function observables as developed 
above) purely on symmetry considerations and, there- 
fore, should be equally applicable to metallic systems. 
Second, a 2D metallic system is "weakly critical" (at dis- 
tances shorter than the localization length) , and the cor- 
responding anomalous dimensions can be studied within 
the perturbativc RG (which is essentially the same as 
Wegner's RG analysis in 2 -I- e dimensions). By a confor- 
mal mapping, a 2D system can be related to the same 
problem in a quasi- ID geometry (with a power-law be- 
havior translating into an exponential decay). Therefore, 
if some symmetry properties of spectra of transport ob- 
servables generically hold at criticality, we may expect to 
see some manifestations of them already in the solution 
of the DMPK equation. 



A. DMPK, localized regime 

In the DMPK approach, the transfer matrix of a quasi- 
lD system is described by "radial" coordinates (w.r.t. a 
Cartan decomposition) Xj, j = 1, 2, . . . , A^, where A^ is 
the number of channels. All transport properties of the 
wire are expressed in terms of these radial coordinates. 
In particular, the dimensionless conductance is 



N 



Tr T = Tr tt^ 



(9.3) 



cosh^ AT, 



(9.4) 



are the transmission eigenvalues, i.e. the eigenvalues of 
T = tt^ (and tH), where t is the transmission matrix. 

The DMPK equations describe the evolution with sys- 
tem length (playing the role of a fictitious time) of the 
joint distribution function for the transmission eigenval- 
ues (or the coordinates Xj), and they have the form of 
diffusion equations on the symmetric space associated 
with the noncompact group of transfer matrices. In the 
localized regime, where the wire length L is much larger 
than the localization length ^, the typical value of each 
transmission eigenvalue becomes exponentially large rel- 
ative to the next one: 1 > Ti > T2 > . . . > Tat. As 
a result, the equations for the random variables Xj de- 
couple, yielding an advection-diffusion equation for each 
Xj. The solution has a Gaussian form, with both the 
average (Xj) and the variance var(Xj) proportional to 
L/^ and with var(A"j) independent of j. Each of the 
symmetry classes therefore gives rise to a set of numbers 
{Xj)/va.j:{Xj) (which depend solely on the corresponding 
symmet ric sp aces) . Remarkably, comparing the above 
results (^)-(|8.1l|) with the known DMPK results, we 



observe that for all symmetry classes one has 



var(A"j 



(9.5) 



(In the case of the chiral classes, we note that Eq. ( |9.5| ) 
holds when the Xj evolve according to the DMPK equa- 
tions with an even number of channels.) 

This result allows us to draw a link between the 
transport quantity Tj and the LDOS observable Vj de- 



get 



4e 



fined in Eq. (5.9). Indeed, if we use the approximation 



which is valid in the localized regime, we 



(T/) - exp{2wg(g + Cj)}, v = var(Xj). (9.6) 

This expression for (T^) has a point q = —Cj/2 of re- 
flection symmetry. [We should add that this requires a 
continuation of Eq. (B^) from its range of actual validity 
to a region of larger q, see the discussion below Eq. (9.2).] 
Now we recall that the scaling of 
the representation (0, . . . , 0, g, 0, . . 
position, see Eq. ( ^.10 ). Hence 



) is determined by 
with q at the j-th 

^ {lyp""), i.e. the 
symmetry point of the multifractal spectrum for Vj is 
exactly —Cj/2. This links Tj with Vj, as stated above. 
We now write 



T — 



T1T2 • • • T„ 



Sn 



T1T2 ■ ■ ■ Tm-1 Srn-1 



, 5™=Ti---T,„, (9.7) 



and draw an analogy between Eq. (9.7) and Eq. (^.g|). 
Specifically, Tm corresponds to Vm (as we have already 
seen earlier) and Sm to Am- To further strengthen the 
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analogy, wc point out that 5*^ can be presented in the 
form of the absolute value squared of a determinant. In- 
deed, eonsider first m = 2. Choose two incoming {p,q) 
and two outgoing (r, s) channels and consider the 2x2 
matrix t^^^ formed by the elements t. 



P,q, J = r,s, 



of the transmission matrix. Then calculate the absolute 
value squared of the determinant of this matrix, and sum 
over the channel indices p, q, r, s: 

^ ^ jdett^^- I — ^ ^ \^pr^qs ips^qr\ 
p.q.r^s p,q,r,s 

= 2(Trttt)2 _ 2Tr(<t^)2 = 2(TrT)2 - 2TrT2 

= 2 [(Ti + T2f - (Tf + r|)] = iTiT2 - 4^2. (9.8) 

The same applies to higher correlation functions: by 
considering the determinant of an m x to matrix t^™^ 
and taking its absolute value squared, we get Sm (up to 
a factor). If the total number of channels is to, this is 
straightforward (the modulus squared of the determinant 
then equals T1T2 • • • T,„); if the total number of channels 
is larger than m, then, strictly speaking, averaging over 
the choice of to channels is required. We expect, however, 
that the determinant will typically behave in the same 
way for any choice of the channels. 

To summarize, the transmission eigenvalues T„i of the 
DMPK model characterize the leading contribution to 
the decay of transport quantities Sm/Sm-i, where the 
Sm are given by the absolute values squared of the de- 
terminants of to X m transmission matrices. There is 
a clear correspondence between the wave-function ob- 
servables vt = Ai/Ai-i and the transport observables 
Tm ~ Sm/Sm-1- In the next subsection we generalize 
this construction to critical systems. 



B. Transport observables at criticality 

We are now ready to formulate a conjecture about the 
scaling of sublcading trans port quantities at criticality. 
It generalizes the relation ( |9.2| ) between the scaling ex- 
ponents of the moments of the conductance (Tq) and of 
the LDOS (xq). 

Consider a system at criticality and take two points 
ri and r2 separated by a (large) distance R. Attach N 
incoming and N outgoing transport channels near each 
of these two points. This yields an N x N transmission 
matrix t. Define B„i as the absolute value squared of the 
determinant of its upper- left m x m corner (i.e., of the 
transmission matrix for the first TO incoming and 
first TO outgoing channels). This lets us build a family of 
transport correlation functions (n < N): 



Mq,q,...qAR} = {Br'^ B^^r'" B!^) 



(9.9) 



where r„ = i?„/i3„_i. The conjecture is that the critical 
index Tg^q^ . g^ determining the leading dependence on R 



of M„ 



wher 



AR) is 



qiq2---qn 



2x, 



qiq2---qn J 



(9.10) 



is the scaling exponent of t he ct- model 



operator (gi, . . . , </„) for the correlator ( 5.10). This is a 
general izatio n of Eq. (3.2). As with Eq. (9.2), the re- 
lation (9.10) is expected to be valid only for qi not too 
large; probably, the condition is qi < —Ci/2 for all i. 



Let us sketch an RG argu ment in favor of Eq. ( 9.10 ) . 
We expect that the quantity ( |9.9D is represented in field- 
theory language as a correlation function of two local op- 
erators (at points ri and r2, respectively), each of which 



has the same scaling properties as i^i 



,q,. 



Performing 



an RG transformation that reduces the scale R down to a 
microscopic scale, we will then get a factor R~'^^iii2 -in . 
After this the correlatio n fun ction becomes of the order of 
unity; thus, we obtain ( 9.10| ). Possibly, a rigorous proof 
may be constructed for class A by a generalization of the 
formula of Ref . . 

It should be stressed that wc do not expect the cor- 
relation functions Eq. (9.£) to show pure scaling: as we 
pointedpaut. not even the moments of the conductance 
show it.E3 



X. CLASSES WITH 0(1) AND U(l) 
ADDITIONAL DEGREES OF FREEDOM 

There are five symmetry classes with cr-model target 
spaces that either have two connected components and 
thus an associated Z2 = 0(1) degree of freedom [classes 
D and Dili], or have Z for their fundamental group due 
to the presence of a U(l) degree of freedom [classes BDI, 
CII, AIII]. These degrees of freedom complicate the ap- 
plication of our Weyl-symmetry argument. 

We mention in passing that the classes at hand are 
the five symmetry classes that feature topological insu- 
lators in ID (precisely because, owing to the 0(1) and 
U(l) degrees of freedom, their cr-model spaces have the 
said topological properties). Below we briefly outline our 
present understanding of the Weyl-symmetry issue for 
these classes and the open questions. 



A. Classes D and Dili 

The target manifolds of the cr-models for these symme- 
try classes consist of two disjoint parts [0(1) = Z2 degree 
of freedom]. In general, the tr-model field can "jump" 
between the two components, thereby creating domain 
walls. The arguments based on the Wcyl symmetry in 
the form presented above apply directly if such domain 
walls are prohibited (i.e. if the cr-model field stays within 
a single component of the manifold). There are several 
situations when this is the case: 

• The DMPK model p£ a quasi-lD wire does not in- 
clude domain walls. E3 This explains the agreement 
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between our symmetry result and the DMPK re- 
sults for these two classes; 



The 0(1) version i 
work model in 2D:I 



the Chalker-Coddington net- 



• A good metal in 2D. (In this case domain walls are, 
strictly speaking, present but their effect is expo- 
nentially small and thus expected to be negligible.) 

Note that the Wcyl-group invariance of the LDOS mo- 
ments for the classes D and Dili yields the symmetry 
point g*/2 = — ci/2 = 0. This implies that the distri- 
bution fu ncti on P{liiiy) is symmetric under Ini/ — > — Ini' 
(see Eq. (1.3) with =0), i.e. Inz^ = is the most prob- 
able (or typical) value. This result is incompatible with 
exponential localization of the eigenstates, which would 
imply exponentially small typical LDOS values. We thus 
arrive at the conclusion that, in the absence of domain 
walls, systems described by the cr-model for class D or 
Dili cannot have a localized phase. The models listed in 
the previous paragraph exemplify this general statement. 

In the case of a good 2D metal in class D or Dili, 
the scaling behavior can be found by perturbative RG, 
with the smallness of the inverse conductance 1/g <^ 1 
ensuring the validity of the loop expansion. In particular, 
the OBerloop RG calculation of the average DOS scaling 
yieldsE3 (i/) cx InX cx g{L). We know that the scaling 
exponents for the LDOS moments depend quadratically 
on q in one-loop approximation (which is governed by 
the quadratic Laplace-Casimir operator). Therefore, in 
view of the q — > —q Weyl symmetry, we expect the LDOS 
moments to behave as 



(10.1) 



It should of course be possible to check this directly by 
a numerical calculation. 



B. Chiral classes 

For the chiral classes, the situation is even more sub- 
tle. We expect that the Weyl-group invariance should 
show up most explicitly in operators that are scalars with 
respect to the additional U(l) degree of freedom. The 
LDOS moments, however, do not belong to this category. 
We leave the SUSY-based classification of operators and 
the investigation of the impact of the Weyl-group invari- 
ance to future work. 



XI. SUMMARY AND OUTLOOK 

In this paper we have developed a classification of com- 
posite operators without spatial derivatives at Anderson- 
transition critical points in disordered systems. These 
operators represent observables describing correlations of 
the local density of states (or wave-function amplitudes) . 
Our classification is motivated by the Iwasawa decompo- 
sition for the (complexification of the) supersymmetric 
CT-model field. The Iwasawa decomposition has the at- 
tractive feature that it gives rise to spherical functions 
which have the form of "plane waves" when expressed in 
terms of the corresponding radial coordinates. Viewed 
as composite operators of the cr-model, these functions 
exhibit pure-power scaling at criticality. Alternatively, 
and in fact more appropriately, the same operators can 
be constructed as highest-weight vectors. 

We further showed that a certain Weyl-group invari- 
ance (due to the Harish-Chandra isomorphism) leads to 
numerous exact symmetry relations among the scaling 
dimensions of the composite operators. Our symmetry 
relations generalize those derived earlier for the multi- 
fractal exponents of the leading operators. 

While we focused on the Wigner-Dyson unitary sym- 
metry class (A) in most of the paper, we have also 
sketched the generalization of our results to some other 
symmetry classes. More precisely, our results are directly 
applicable to five (out of the ten) symmetry classes: the 
three Wigner-Dyson classes (A, AI, All) and two of the 
Bogoliubov-de Gennes classes (C and CI). Moreover, they 
should also be valid for the remaining two Bogoliubov-de 
Gennes classes (D and Dili), as long as cr-model domain 
walls are suppressed (i.e. the cr-model field stays within 
a single component of the manifold). Our results imply 
that in this situation the system is protected from An- 
derson localization. In other words, localization in the 
symmetry classes D and Dili may take place only due to 
the appearance of domain walls. 

We have further explored the relation of our results 
for the LDOS (or wave-function) correlators to transport 
characteristics. We have constructed transport observ- 
ables that are counterparts of the composite operators 
for wave-function correlators and conjectured a relation 
between the scaling exponents. 

Our work opens a number of further research direc- 
tions; here we list some of them. 

(i) Verification of our predictions by numerical simu- 
lation of systems housing critical points of various 
dimensionalities, symmetries, and topologies would 
be highly desirable. While the LDOS multifractal 
spectra have been studied for a considerable num- 
ber of critical points, the numerical investigation of 
the scaling of subleading operators is still in its in- 
fancy. Preliminary numerical results for the spectra 
of scaling exponents of the momcHte (Aj) and (A|) 
at the quantum Hall critical pointuJ do support our 
predictions. Furthermore, it would be very inter- 
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esting to check numerically our predictions for the 
scaling of transport observables. 



Pi P2 P3 Pa Pb P6 



Pi P2 P3 Pi P5 P6 



(ii) As mentioned in Sec. X B , it remains to be seen to 
what extent our results can be generalized to the 
chiral symmetry classes, and what their implica- 
tions for observables will be. 

(iii) In this work, we have studied critical points of non- 
interacting fermions. In some cases the electron- 
electron interaction is RG-irrelevant at the fixed 
point in question, so that the classification remains 
valid in the presence of the interaction. An example 
of such a situation is provided by the integer quan- 
tum Hall critical point fiitb. a short-range electron- 
electron interaction.OO However, if the interac- 
tion is of long-range (Coulomb) character, the sys- 
tem is driven to another fixed point. (This also hap- 
pens in the presence of short-range interactions for 
fixed points with spin-rotation symmetry: in this 
case, the Hartree-Fock cancelation of the leading 
term in the two-point function (1.5) docs not take 
place.) The classification of operators and relevant 
observables at such interacting fixed points, as well 
as the analysis of possible implications of the Weyl- 
group invariance, remain challenging problems for 
future research. 
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Appendix A: Young diagrams, tableaux, and 
symmetrizers 
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FIG. 2. Young diagram A = (6^ 4^ 2, 1) = [6,5,4^,22] (left) 
and its conjugate A = (6, 5, 4^ 2^) = [6^, 4^, 2, 1] (right). 



agram A = (6,6,4,4,2,1) 
4 distinct row lengths. 



^ 4^, 2, 1) with 6 rows and 



At the top of the diagram in Fig. || we also display 
the numbers of boxes pi,p2, ■ ■ ■ ,Pm in each column. Like 
the numbers , . . . , q„ , these completely specify the di- 
agram, and we will use (for the same diagram) the al- 
ternative notation A = [pi,p2, ■ ■ ■ ,Pm] = [^iS ■ • ■ > ^s'Ij 
where the second notation means that the partition by 
the integers pi has aj copies of the integer kj, 1 < j < s. 
In this notation the diagram shown in Fig. on the left 
is A = [6,5,4^,2'^]. The numbers pi also define the con- 
jugate diagram A = (pi,p2, • ■ • ,Pm)- For illustration, the 
diagram A = (6, 5, 4^, 2^) = [6^, 4^, 2, 1] is shown in Fig. | 
on the right. (Notice that the number s of distinct parts 
is the same for a Young diagram and its conjugate, and 
is the same as the number of "corners" on the boundary 
of the diagram.) The number of boxes of A and A, called 
the size of A, is denoted by 



N. 



(Al) 



For a given Young diagram, the integers h and are 
related by 



ai 
a2 



''Si 
ls-1 



In this paper we use a standard notation for Young 
diagrams, see, for example, Ref. [Q. Thus, the Young 
diagram corresponding to the partition 91 + 92 + ■ . ■ + 9n 
(where the integers qj are subject to qi > q2 > ... > 
9n > 0) is denoted by A = {qi, 52, ■ . • , 9?i), and consists of 
left-aligned rows with the top row containing qi boxes, 
the next row containing q2 boxes, etc. Another notation 
that we will use is A = (Zj\ . . . , i^') to denote the parti- 
tion that has bi copies of the integer h, I < i < s. For 
an example, on the left in Fig. ^ we show the Young di- 



flj — Zi — ^2- 

Solving this for gives 

h = ai + a2 + . . . + Us, 
h — ai + a2 + ■ ■ ■ + tts-i, 

Is = ai. 



(A2) 



(A3) 
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Similar relations exist between h and bf. 



bi 
b2 

bs 

h 

k2 



ks—i k 

ki - k2. 

&i + 62 4 
&i + &2 4 



(A4) 



+ bs, 
+ bs-i, 



kg = bi. 



(A5) 



Young diagrams are used to label irreducible repre- 
sentations (irreps) of the permutation groups and some 
classical matrix groups. Irreps of G = U(n) of polyno- 
mial type are in one-to-one correspondence with Young 
diagrams that have at most n rows. The eigenvalue of 
the quadratic Casimir operator in the irrcp of U(n) with 
Young diagram A = (gi, 92, ■ • ■ , is 



a2 



(A,n) = ^qi{q^ + Jl+l-2i). 



(A6) 



It is known that the quadratic Casimir eigenvalue for the 
conjugate Young diagram A is related to the one for A by 



a2(A,?i) = -a2(A, 



(A7) 



Next we need the notion of Young tableaux. A tableau 
T is a Young diagram A with each of its boxes filled with 
a positive integer from the set {l,2,...,iV = |A|}. A 
tableau is called semistandard if the integers in the boxes 

i) weakly increase from left to right along each row and 

ii) strictly increase from top to bottom along each col- 
umn. The minimal semistandard tableau (denoted by 
Tmin) for a given shape A = [fc"^ , . . . , fc""] is the one where 
all the integers in the first row are 1, in the second row 
2, and so on, up to ki in the last row. A semistandard 
tableau is called standard if it is filled according to the 
above rules so that each number from the set {1, . . . , N} 
occurs exactly once. A normal tableau (which we denote 
by To) is a standard Young tableau in which the numbers 
arc in order, left to right and top to bottom. If a tableau 
T is obtained by filling a Young diagram A, we say that A 
is the shape of T. To give an example. Figure ^ shows a 
semistandard, the minimal, a standard, and the normal 
tableau of shape (6, 4^, 2). 

The permutation group Sn acts on tableaux with N 
boxes by permuting the integers in the boxes. If cr e Sn, 
we denote by aT the tableau which has the number a{i) 
in the box where T has i. For each Young diagram A of 
size N we define i?(A) and C(A) as the subgroups of Sn 
that preserve the rows and columns of A, respectively. 
One can consider formal linear combinations of the el- 
ements of Sn (these form what is known as the group 
algebra A{Sn) of Sn) and define the following opera- 
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(a) Semistandard 
tableau 



(b) Minimal tableau 
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(c) standard tableau 



(d) Normal tableau 



FIG. 3. Young tableaux 



tors: 

aA= 51 ^' H sgn(r)T. (A8) 

(T6i?(A) t6C(A) 

When acting on a tableau T, the operator a\ symmetrizes 
all the numbers in T along its rows. Similarly, b\ anti- 
symmetrizes entries of a tableau along the columns. Fi- 
nally, the Young symmetrizers are defined as the prod- 
ucts 



(A9) 



Sometimes one uses an alternative definition of the Young 
symmetrizers where the order of the operations of sym- 
metrization and antisymmetrization is reversed: 



CA flA^A- 



(AlO) 



All operators oa, b\, c\, and ca are idempotent; this 
means that their squares are proportional to the opera- 
tors themselves: 



nRa\, 



b\ = ncb\, 



(All) 



where and nc are the orders of the subgroups i?(A) 
and C(A); and n\ is another positive integer. 

For Young diagrams of type (g) = the Young sym- 
mctrizer reduces to the total symmetrizer along the single 
row. Similarly, for the diagrams of type (P) = [p\ the 
Young symmetrizer is the total antisymmetrizer along 
the single column: 



=(9) 



C[p] 



(A12) 



To illustrate these operators, consider the normal 
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tableau Tq for A = (2, 1). In that case, 
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(A13) 
(A14) 

(A15) 

(A16) 



For the purposes of this paper we also need to con- 
sider tableaux filled by points and wave-function sym- 
bols rather than integers. The action of the permutation 
groups and Young symmetrizers on such tableaux is de- 
fined in the same way as on tableaux filled by integers: 
the points and wave functions are permuted according 
to their positions in a tableau. Now, suppose we have a 
young diagram A and two tableaux, and of shape 
A, filled with wave functions and points, respectively. We 
can define a pairing of these two tableaux as the following 
product of wave functions: 



(A17) 



where i runs over the boxes of the diagram A. The 
tableaux used in this definition need not be standard or 
semistandard but can be arbitrary. For example, for the 
following two tableaux 



1p2 












1p3 
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r2 




ri 






r4 







(A18) 



the corresponding product of wave functions is 

*(3,2,i)(Tv„T,) ^ V^i(ri)V'2(ri)V3(r4)V'4(r2)V5(r3)- 



When an element s £ ^(5*1^1 ) of the group algebra of 
the symmetric group S\\\ acts on one of the arguments 
of '^x[T^^Tr)^ we understand, say 



(A19) 



as the linear combination of the corresponding products. 
It is easy to derive some useful properties of such actions. 
First of all, it is clear that if we permute the entries in 
both tableaux and in the same way, then we do 



not change the pairing of these two tableaux: 

«'A(aT0,tTT,) = *a(T^0,T,), a £ S\x\- (A20) 

Now take the inverse cr^^ of a permutation a e 'S'|a|, and 
apply it to both tableaux in ^a(2^0, crTr). Then 



(A21) 



Next, if the permutation a in the last equation runs over 
either of the subgroups i?(A) or C{X), its inverse does 
the same. Moreover, the paritie s of a and are the 
same. Therefore, summing Eq. (A21) over R{X) or over 
C(A) with appropriate sign factors, we get 

■^x{T^,axTr) = -^x{axT^,Tr), (A22) 
■9xiT4„bxTr)^-^xibxT.^,Tr). (A23) 

Finally, using the last two equations it is easy to obtain 



^xiT^,cxTr) = ^xibxT^,axTr) = «'a(5aT^, T,), 



(A24) 

^x{T^„cxTr) = ^x{axT4.,bxTr) = 1'a(caT^, T,). 

(A25) 

The combinations of wave functions that play a special 
role in the paper are ^a(T, caT), where T is a standard 
tableau of shape A. It is these particular combinations 
that lead to pure scaling operators in the cr-model, see 
Section VI. For example, we can take both and Tr to 



be the normal tableau and obtain ^'a(To, caTq). 

Let us look at a few simple examples. If we take both 
and Tr to be the normal tableau for the diagram (2,1) 
and act on Tr by th e Young symmetrizer C(2.i) of ( |AS| ), 
then by using (A15) we get 



*(2,i)(To,C(2,i)To) = ■i/'i(i"i)V'2(r2)V'3(r3) 

- V'i(r3)^2(r2)V'3(ri) + '0i(r2)-!/'2(ri)V'3(r3) 

- Vi(r3)^/'2(ri)^/'3(r2) 

= i^2(r2)i?2(ri,r3) + i^2(ri)i:'2(r2, rg). (A26) 

It is an easy exercise to show that one has the same ex- 
pression for ^'(2,i)(c(2,i)To,To). 

In the notation of Eq. ( A19| ) the Slater determinants 
( 4.16 ) can be written as 



Dp(ri, . 



^(ip)(?o, ^(ip)To)- 



(A27) 



If we build our product of wave functions by taking 
the minimal semistandard tableau of a given shape A = 
[fc"\ . . . , kg"] for both T^, and Tr, then the operation of 
symmetrization along the rows is clearly redundant, and 
we get 



*A(rinin,CAT'i„in) OC x{Tmin, bxT^^iri) 

= *A(&AT„,i„, r„,i„) = Dii ■ ■ ■ Dii , 



(A28) 



where each determinant Dj is evaluated on wave func- 
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tions ipi, . . . at the points ri , . . . , . Adopting the 
notation A = (51, 52, ■ • • , <Zn), an alternative form of this 
expression is 



(A29) 



Appendix B: Construction of highest-weight vectors 

In th is ap pendix we construct the scahng operators 
(3.21), ( |6.25 ) by using the idea of the highest- weight vec- 
tor sketched in Sec. VI F. As was discu ssed there, we first 
focus on hnear functions fXyiQ) ( |6.35| ) of the matrix el- 
ements of the cr-model field Q specified by a matrix Y: 
I^y{Q) = Tr(y(5). Let us denote by Eij the matrix which 
contains the number one at the intersection of the i-th 
row with the j-th column and zeros everywhere else. Such 
matrices are sometimes called "matrix units" . Individual 
matrix elements of Q can be written as 



/iij., (Q) = Q] 



(Bl) 



We begin with two simple examples: functions on a 
sphere S"^ and on a hyperboloid ; these arc symmetric 
spaces of compact and non-compact type, respectively. 



1. Functions on a sphere 

Consider the space of functions on the two-sphere 

G/K = U(2)/U(l) X U(l) = SU(2)/U(1) = S^. (B2) 

To make the presentation here similar to the general case 
considered later, we represent points on the sphere by the 
matrix 



= .9A.g \ 



(B3) 



where g £ SU(2) and A = (73 is the third Pauli matrix. 
Using a parametrization of SU(2) by Euler angles. 



g = 
we get 



-i(0+VO/2 cos ^ „e-*('^-'^)/2 f 
gz(0-V.)/2 sin I e'(0+VO/2 , 



Q 



X3 Xi - 1X2 

Xi+iX2 -X2, , 



(B4) 



(B5) 



the Cartan generator is X3, and the raising and lowering 
generators (in the complexiflcation sl(2,C) of su(2)) are 
X± = iXi T X2. Notice that 



A+ — £^12, A_ — E21- 



Now the function 



iQ) = Mx+ (Q) = Q21 = xi+ 1X2 = sin ( 



(B7) 



(B8) 



is a hig hest-w eight vector for the SU(2) action. Indeed, 
by Eq. (S.36) we have 



X+'fi = fJ-[x+.x+] = 0, 
Xaipi = ^J.[X3.x+] = y'-ix+ = -iv>i- 
Powers of this function, 



(B9) 
(BIO) 

(Bll) 



are also highest- weight vectors. To make them globally 
well defined on the sphere, the power / has to be a non- 
negative integer. Constant multiples of ipi arc known as 
the spherical harmonics Yu in the irreducible representa- 
tion of SU(2) of dimension 2Z -I- 1. 

There exist other choices of Cartan and nilpotent sub- 
algcbras. For example, if we chose Xi as the Cartan 
generator, and X'_^ = —A3 -I- iX2 as the raising operator, 
then the (linear) highest-weight vector would be 



^'1 



^{Qll-Q21+Ql2-Q22) 

i{x3 + 1x2) ~ icos9 — sin 6* sin 0. (B12) 



Similarly, the choice of A2 as the Cartan generator would 
lead to the highest-weight vector 



ip'( ~ X3 + ixi = cos 9 + i sin cos (f>. 



(B13) 



Both if'i and can be raised to non-negative integer 
powers to produce other highest-weight vectors. 

In this example of functions on a compact symmetric 
space, all three choices of Cartan subalgebra or highest- 
weight vector are equivalent and can be transformed into 
each other by an element of G = SU(2). (In fact they are 
just three "points" on an SU(2)-orbit of Cartan subalge- 
bras or highest-weight vectors.) This will not be the case 
in our next example of functions on the two-hyperboloid. 



where 



2. Functions on a hyperboloid 



xi — sinO coscj), X2 = sin 6* sin 0, X3=cos6' (B6) 

are the three basic functions which arise by restricting 
the Cartesian coordinates of the Euclidean space M.^ to 
the sphere x\ + x\+ x^~l. 

Let us choose Xk = ak/{'2-i) as our system of gener- 
ators of the Lie algebra su(2). The standard choice for 



We now consider the space of functions on a non- 
compact analog of the two-sphere, the two-hyperboloid 

G/A' = SU(1,1)/U(1) = i72. (B14) 

One may view this space as a non-compact variant of 
the sphere S"^, by analytically continuing the compact 
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angle 9 to the non-compact radial variable on (de- 
noted by the same symbol 9). If we make the replace- 
ment 9 ^ iO in the function (|B8[), we get isinh6'e*'^. 
While this function is a highest-weight vector for some 
choice of the Cartan subalgebra, it is not positive on 
the hyperboloid, so it cannot be raised to an arbitrary 
complex power. However, there exist other, inequivalent 
choices of the Cartan subalgebra which do give the de- 
sired positivity property. In f act, if we analytically con- 
tinue 9 i9 in, say, Eq. ( B13 ), we get the highest- weight 
vector cosh 9 — sinh 9 cos </>, which is strictly positive on 
and, therefore, can be raised to an arbitrary complex 
power. Here is how it is done more formally. 



Matrices g G SU(1, 1) satisfy the relation 



(B15) 



and can be parametrized in terms of generalized Euler 
angles as 



e'('^+'^)/2 cosh f -ie*W-^)/2sinhf 
•g-^(0-^)/2 gij^h I e-'('^+'^)/2 cosh § 



(B16) 



Elements of the coset space G/K = are represented 
by matrices 



Q = 9<^39 ^ 



X3 
ixi - X2 



(B17) 



The matrix elements xx^x^TX-i may be viewed as the 
Cartesian coordinates of the Euclidean space re- 



stricted to the hyperboloid x\ - 



the parametrization (B16) we express them as 



1. By adopting 



x\ = sinh cos X2 — sinh sin 0, x'3 = cosh^ 



(B18) 



The Lie algebra su(l, 1) ~ is spanned by the ma- 
trices iX\^ 1X2, and X3. Choosing for the Cartan 
generator and X± = iXi =F X2 for the nilpotent genera- 
tors, we get the highest- weight vector 



Vi ~ Q21 = ixi — .T2 = i sinh( 



(B19) 



This is the analog of (B8) for the hyperboloid, and it 
is not a positive function. To obtain a positive highest- 
weight vector, we need to choose a linear combination of 
iXi and 1X2 for the Cartan generator. 



Thus let the Cartan generator be iXi cos a — 1X2 sin a 
for some choice of parameter a. Taking 

X'j^ =^ X2 cos a + Xi sin a + iX^ 

for the raising operator, we have the following expression 



for the corresponding highest- weight vector: 



Mx; = \ {Qii - (522 + e ' 
x^ + X2 cos a + xi sin a 
cosh 6* + sinh6'sin((/) + a), 



{21) 



(B20) 



which already arose in the closely related context of Eq. 
(4.11). An arbitrary complex power of this positive func- 



tion is also a highest- weight vector; 
(/3^ = (cosh^ + sinh6'sin((?!) + a))', 



(B21) 



and this function (or rather, its extension to the SUSY 
setting) is the cr-model scaling operator for the g-th power 
of the local density of states. 

While it is clear by inspection that the function Lp'^ is 
positive, a more formal proof that generalizes to cases of 
higher rank is as follows. We write 

ip'i — iTr[(cr3 — i<72 cos a — idi sin a) gcrzg"^] 

= iTr[(l - (Ti cosa + (73 sina).g.g'^], (B22) 



where we have used the SU(1, 1) defining relation (B15). 
The matrix H = (1 — cti cos a + a2 sin a) /2 is a projection 
operator: H^' = H = H^. We thus see that the function 
Lp'^ is the manifestly positive expectation value of gg^ > 
in the eigenvector of the projector H with eigenvalue 1. 



3. Arbitrary n, compact case 

We now come back to the general case of functions on 
the compact symmetric space for class A: 



G/K = U(27i)/U(n) X U(n), 



(B23) 



which arises from the use of fermionic replicas. Elements 
of this coset space, or points on the manifold, are repre- 
sented by matrices Q = gKg~^, where now 



A 




(B24) 



We begin with a choi ce of root-space decomposition 
= n+® f)©n_, see Eq. ( 6.30| ). We take \) to be spanned 
by the diagonal matrices and n+ (ri-) be spanned by the 
upper (respectively, lower) triangular matrices in g = 
gl(2??., C). Schematically, 



/* 

* 





Vo 



o\ 





* 
*/ 









Vo 



* 
0/ 



(B25) 



We also need the (refined) Cartan decomposition g 
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p+©fi©p_, where t = gl(n, C) © gl(n, C) is the complex- 
ified Lie algebra oi K = U(n) x U(n), while p± are the 
eigenspaces of the adjoint (or commutator) action of S3: 



* 



* 
* 



p+ 



\0 ... * ... */ 

Note the commutation relations 




Vo ... ... 0/ 

(B26) 



[P+,P+] =0, 



(B27) 



as well as the following decomposition of the space n+ of 
raising operators: 



n+ =P+©(n+n«). 



(B28) 



Our attention now focuses on t he space of complex- 
valued functions /iy (sec Eq. ( 3.35| )) for Y G p+. We will 
use such functions as building blocks to construct func- 
tions that have the properties of a highest-weight vector, 
see ( |6.31 ). In fa ct, by the second set of commutation re- 
lations in (B27) any function /iy for F G p+ is already 
annihilated by all first-order differential operators that 
represent generators of p+: 



= Ai[x,y] =0 ioT X,Y e p^ 



(B29) 



However, fiy for general y € p.). is not annihilated by all 
raising operators from n+ H 5. To implement this annihi- 
lation condition, we construct certain polynomials of the 
matrix elements of Q in the following way. 

For 1 < i, j < n we introduce the functions 



Notice that the matrix E, 



Q 



2n+l-j,i- 



(B30) 



'i,2n+i~j G P+, SO the functions 
Vij are exactly of the type discussed in the previous para- 
graph. Now we will demonstrate that for any integer m 
in the range 1 < m < n the m x m determinant 



Dot 



I'll 



(B31) 



is a highest- weight vector for the decomposition (5.30). 

We first establish that Xf„i = for all X G n+. Due to 
the decomposition (|B28| ) there are two cases to consider. 
First, let X G p+. Then, as we have already mentioned, 
Xvij — and, therefore, X/,„ = 0. Now let X be in the 
space n+ n 6, which is spanned by Eai ioi 1 < i < i' < n 
and Ejji for n -|- 1 < j < j' < 2n. If X — Eui with 



m < i' < n, then we still have Xfm = 
1 < < m the commutator 



0, since for all 



(B32) 



vanishes due to i' 7^ i" . Now let X ~ Eui with 1 < i < 
i' < m. In this case we obtain 



dt 
d_ 

dt 



Dot [Tr( 



'E,„ 



(B33) 



The matrix under the determinant sign in the last equa- 
tion factorizes as (1 -I- tEi'i)v, therefore 



Eii' fni — fjr 



dt 



t=o 



Dct{l + tE,,,) = 0, 



(B34) 



since the determinant in the last equation does not de- 
pend on t. 

It remains to show that /,„ is an eigenfunction of the 
operators from f). To this end we express H G t) as H = 
Y^iiZi^iEii, where hi G C. This is a diagonal matrix, 
and so is e*^ = X]i=i Then we have 



Huij = — I > e 



''^Tv{EkkE,^2n+i-jEuQ) 



dt\t=o 

and, indeed, the property ( |6.31 ) follows: 
d 



(B35) 



- Jt 



dt \t=Q 



(B36) 



Since all functions /,„ for 1 < m < n have the highest- 
weight property, so does the product 

_ f?i-<?2 fq2-q3 fqn-i~qn fq„ mQ7A 
<<5(gi,...,g„) - 7l h ■■■Jn-l Jn (^-^ ' ) 

for a weakly decreasing sequence of n integers gi > (72 > 
■ ■ ■> <ln > 0. The powers in this expression are restricted 
to be no n- negative integers, since the functions /,„ are 
complex-valued. The functions i^(gi,...,g„) are the most 
general highest-weight vectors in the present situation. 

Note that if iJ = T^^^i hiEu is a diagonal generator 
of GL(n,C), then 



i=l 



)<^(ql,...,9„)• 



(B38) 



By standard facts of representation theory it follows that 
(gi , . . . , g„) my be interpreted as the sequence of numbers 
determining the Young diagram of an irreducible repre- 
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sentation of GL(n, ( 



triangular. As before, we introduce a set of functions 



4. Arbitrary n, non-compact case 



We now turn to the general non-compact situation for 
class A and consider the space of functions on 



G/K = \]{n,n)/\]{n) x U(ri), 



(B39) 



which results from the use of bosonic replicas. Points on 
this manifold are represented by matrices Q = gAg^^, 
where A = S3 is defined in Eq. (B24). Elements of the 
pseudo- unitary group U(n, n) satisfy 
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(B40) 



Hence the functions fJ-viQ) can be rewritten as 

^iY{Q) = TiiYg^^g-') = TT{gg^j:3Y). (B41) 

We will try to follow the development of the compact 
case as much as possible. One major change comes from 
the fact that by the hyperbolic nature of the Lie algebra 
= u{n,n), there exist several G-inequivalent choices of 
Cartan subalgebra t). For our purposes, the good choice 
to consider is as follows. 

We make an orthogonal transformation of the standard 
basis {ci} of C^" to introduce a new basis {ii}: 



Cj + e 



j+n 



V2 



e2n+l-j 



V2 



(B42) 



where 1 < j < n. We then define linear operators Eij 
on C^" by the relation Eijik = Sjkei. Thus, in the new 
basis the operators Eij are the "matrix units" . They can 
be expressed in terms of the matrix units with respect to 
the original basis: 

E2n+l-i.j = ^{^ij ^ Ei^n.j + Ei^j^n ^ -E'.i+,i,j+„ ) , 
Ei^2n+l-j = \{Eij + Ei^n.j ~ Eij^n ^ Ei-^n,j+n), 
^2ri+l-i,2ri+l-j = ^{Eij — Ei^n,j ^ Ejj-^n + Ei^nJ+n)- 

With these conventions, let us choose the Cartan sub- 
algebra f) and the subalgebra n+ of raising operators as 
follows: 

I) = spanc {E'jj , 1 < i < 2n} , (B43) 
n+ = spaiicj-Ey , 1 < « < j < 2n}. (B44) 

Thus our Cartan generators in the transformed basis e 
are still diagonal and the raising operators are still upper 



= Q 



z+n Qj-\-n,i Q j-^7i,i-\-n ) (B45) 

for 1 < i,j < n. We point out that in the advanced- 
retarded space this is exactly the structure that appeared 
before in Eq. 

We now define functions /,„ for 1 < m < n in the same 
way as in the compact case: 



(m) 



I'll 



(B46) 



By the same argument as for the compact case, each 
of the fm has the properties (6.31) of a highest-weight 
vector. 

Moreover, each of the functions fm is real- valued and 
positive. This is s een as follows. Recalling the second 
expression in ( B41 ) we have 



Tr(.g<7^n, 



where 



'i^d,Ei^2n+l-j 



^2n+l-i,2n+l-j- 



(B47) 



(B48) 



This shows that uij is, in fact, a single matrix element of 
the positive definite matrix gg'^ in the new basis, in the 
subspace spanned by {e2,i+i-j}i<i<„: 



55' 



2n+l-j,2n+l-i ' 



(B49) 



Then the determinant = Det[i^^"'^] is a principal mi- 
nor of this positive definite matrix, and, therefore is pos- 
itive as well. Hence, in the construction of a general 
highest- weight vector. 



n 



qi-q2 ^92 -93 
J2 



' ' n—1 J n 



(B50) 



we may take the qi to be arbitrary complex numbers. 

Notice, on the other hand, that the functions fm are 
the principal mi nors of the appropriate block of the ma- 
trix Q, see Eg. (p345| ). Therefore, in the notation intro- 
duced in Sec. |VI| for these minors, /,„ = c?,„, the highest- 
weight vecto r is the same as the function tpq^a 
in Eq. 

Finally, we comment that the generalization to the su- 
persymmetric case is straightforward. We simply need to 
replace all traces by supertraces and determinants by su- 
perdeterminants. Otherwise, everything goes through in 
the same way as before. For the purposes of this paper, 
it is sufficient to consider only the non-compact (boson- 
boson) sector of the super cr-model. In this case the pow- 
ers qi can again take arbitrary complex values. We thus 
reproduce eigcnfunctions ( |6.21 ) which, as explained in 
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Sec. VI D, are none other than the A''-radial functions 
( 3.25 ) of the Iwasawa-decomposition approach. If we do 
not restrict ourselves to the boson-boson sector, we ob- 
tain a broader class of eigenfunctions that, b y the same 
token, will be equivalent to the plane waves (3.15). The 
powers pi corresponding to the compact sector are then 



non-negative integers as in Appendix B 3 
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TABLE II. cr-model spaces. The cr-model target spaces for the localization problem fall into the large families of Riemannian 
symmetric superspaces. The last two columns list the compact and non-compact components of their underlying manifolds. 



bymmetry 
Class 


iNLuJVl 
(n-c|c) 


Compact space 
(// sector) 


Non-compact space 
(bb sector) 


A (UE) 


AIII 1 AIII 


U(2n)/U(n) x U(n) 


U(n,n)/U(n) x U(n) 


AI (OE) 


BDI|CII 


Sp(4n)/Sp(2n) x Sp(2n) 


0(n,n)/0(n) x 0(n) 


All (SE) 


CII|BDI 


0(2n)/0(n) x 0(n) 


Sp(2n,2n)/Sp(2n) x Sp(2n) 


AIII (chUE) 


A|A 


U(n) 


GL(n,C)/U(n) 


BDI (chOE) 


AI|AII 


U(2n)/Sp(2n) 


GL(n,]R)/0(n) 


CII (chSE) 


AII|AI 


U(n)/0(n) 


GL(n,H)/Sp(2n) 
= U*(2n)/Sp(2n) 


C (SC) 


DIII|CI 


Sp(2n)/U(n) 


S0*(2n)/U(n) 


CI (SC) 


D|C 


Sp(2n) 


SO(n,C)/SO(n) 


BD (SC) 


CI|DIII 


0(2n)/U(n) 


Sp(2n,M)/U(n) 


Din (SC) 


r|D 


0(n) 


Sp(2f/.C)/Sp(2/7) 



TABLE III. Root systems for the cr-model target spaces. We choose the system of positive roots such that, in the notation used 
in the table, 1 < j < k < p and 1 < m < I < r (notice the opposite choice for the bb and // sectors). As appropriate for our 
a-modol target spaces, wo only consider the orthogonal groups in even dimensions. The last two columns list the coefficients 
of the expansions of the half-sum of positive roots p = X)^=i ^j^j + * X)i=i ^'2/'- 



Symmetry 
Class 


NLctM 
(n-c|c) 


Xj Xk 


Xj ~\~ Xk 


2xj 


iiUl - Urn) 


i{Vl + Vm) 


2iy; 


Xj - iyi 


Xj + iyi 


P 


r 




bi 


A 


AIII 1 AIII 


2 


2 


1 


2 


2 


1 


-2 


-2 


N 


N 


l-2j 


21 - 1 


AI 


BDIjCII 


1 


1 





4 


4 


3 


-2 


-2 


2N 


N 


-.7 


AI - 1 


All 


CII|BDI 


4 


4 


3 


1 


1 





-2 


-2 


N 


2N 


3 - 4 j 


/ - 1 


AIII 


A|A 


2 








2 








-2 





N 


N 


l-2j 


21- 1 


BDI 


AI|AII 


1 








4 








-2 





2N 


N 


I 


4/- 2 


CII 


AII|AI 


4 








1 








-2 





N 


2N 


2-Aj 


'4 


C 


Dili 1 CI 


4 


4 


1 


1 


1 


1 


-2 


-2 


N 


2N 


l-Aj 




CI 


D|C 


2 


2 





2 


2 


2 


-2 


-2 


N 


N 


-2,? 


2/ 


BD 


CI 1 Dili 


1 


1 


1 


4 


4 


1 


-2 


-2 


2N 


N 


1-J 


4/ - 3 


Dili 


C|D 


2 


2 


2 


2 


2 





-2 


-2 


N 


N 


2-2i 


2/ -2 



